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A  method  for  treating  nonlinear  stochastic  systems  is  described  which 
it  is  hoped  will  be  useful  in  both  the  quantum-mechanical  many-body  problem 
and  the  theory  of  turbulence.   In  this  method  the  true  problem  is  replaced 
by  models  that  lead  to  closed  equations  for  correlation  functions  and  aver- 
aged Green's  functions.  The  model  solutions  are  exact  descriptions  of  pos- 
sible dynamical  systems,  and  as  a  result  they  display  certain  consistency 
properties.  For  example,  spectral  components  of  Green's  fiinctions  which  must 
be  positive-definite  in  the  true  problem  automatically  are  so  for  the  models. 
The  models  involve  a  new  stochastic  element:  Random  couplings  are  introduced 
among  an  infinite  collection  of  similar  systems,  the  true  problem  corresponding 
to  the  limit  where  these  couplings  vanish.  The  method  is  first  applied  to  a 
linear  oscillator  with  random  frequency  parameter.  The  mean  impulse-response 
function  of  the  oscillator  is  obtained  explicitly  for  two  successive  models. 
The  results  suggest  the  existence  of  a  sequence  of  model  solutions  which  con- 
verges rapidly  to  the  exact  solution  of  the  true  problem.  Applications  then  are 
made  to  the  Schrodinger  equation  of  a  particle  in  a  random  potential  and  to 
Burgers'  analog  for  turbulence  dynamics.  For  both  problems,  closed  model  equa- 
tions are  obtained  which  determine  the  average  Green's  function,  the  amplitude 
of  the  mean  field,  and  the  covariance  of  the  fluctuating  field.  The  model  solu- 
tions can  be  expressed  as  sxims  of  infinite  classes  of  terms  from  the  formal 
perturbation  expansions  of  the  solutions  to  the  true  problems.  It  is  suggested 
that  correspondence  to  stochastic  models  may  be  a  useful  criterion  to  help 
judge  the  validity  of  partial  summations  of  perturbation  series. 
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1.     INTRODUCTION 

This  paper  is  intended  to  introduce  a  method  for  treating  certain  problems 
where  the  dynamical  equations  are  nonlinear  in  stochastic  quantities.     The 
quantum-mechanical  many-body  problem  *     and  the  theory  of  turbulence  *     are  two 
fields  of  current  interest  where  it  is  hoped  that  the  method  will  prove  useful. 
In  such  problems,  there  arise  from  the  dynamical  equations  an  infinite  hierarchy 
of  coupled  equations  which  relate  given  ensemble  averages  to  successively  more 
complicated  ones.     An  equivalent  statement  is  that  the  prediction  of  a  given 
average  over  a  finite  time  requires  the  initial  knowledge  of  an  infinite  number 
of  averages.     This  situation,   which  commonly  is  called  the  closure  problem, 
arises  even  when  the  nonlinear  stochastic  terras  are  linear  in  the  dynamic  vari- 
ables.    An  example  is  linear  wave  propagation  in  a  medium  with  random  refractive 
index  fluctuations.       Here  the  equation  for  the  ensemble-averaged  wave  amplitude 
forms  the  base  of  an  hierarchy  involving  successively  higher  cross-moments  of 
the  joint  distribution  of  index  and  amplitude  fluctuations. 

A  formal  solution  to  the  (fynamical  equations  of  any  of  the  problems  men- 
tioned above  may  be  obtained  by  treating  the  nonlinear  terms  as  a  perturbation 
and  expanding  by  iteration.   *       One  may  then  approximate  statistical  quantities 
by  either  truncating  this  expansion  or  summing  tractable  classes  of  terms  to  all 
orders.     Another  (and  related)  approach  is  to  discso'd  all  cumulants  of  the  sta- 
tistical distribution  above  a  certain  order.     Then  all  averages  are  expressible 

in  terms  of  averages  of  this  order  and  below,   thereby  providing  a  closure  of  the 

on 
hierarchy  of  coupled  statistical  equations.   ' 

In  the  method  to  be  presented  here,  the  true  problem  is  replaced  by  models 

that  lead, without  approximation,  to  closed  equations  for  correlation  functions 

and  averaged  Green's  functions.     The  model  solutions  are  exact  descriptions  of 


possible  dynamical  systons,  and  consequently  they  have  certain  consistency 
properties  which  can  be  lacking  in  the  approximation  schemes  mentioned  above. 
For  example,   spectral  components  of  Green's  functions  which  must  be  positive- 
definite  in  the  true  problem  automatically  are  so  in  the  models.     A  related 
property  is  that  covariances  satisfy  certain  realizability  inequalities. 

The  models  are  constructed  by  introducing  dynamical  couplings  among  an 
infinite  collection  of  similar  systems,     the  true  problem  corresponding  to  the 
limit  in  which  these  couplings  vanish.     The  coupling  coefficients  change  ran- 
domly from  one  individual  system  in  the  collection  to  another.     Thus  they  con- 
stitute a  new  stochastic  elanent  not  present  in  the  true  problem.     The  models 
are  most  easily  formulated  in  terms  of  a  collective  representation  in  which  the 
variables  are  linear  combinations  of  those  of  all  the  individual  systems. 

The  closed  statistical  equations  which  characterize  the  models  are  obtained 
by  averaging  over  an  ensemble  of  realizations  of  the  collection  of  coupled  sys- 
tems.    When  iteration  expansions  are  generated  for  the  averages  of  basic  i,nterest, 
it  is  found,   using  the  collective  representation,   that  the  random  couplings 
result  in  the  cancellation  of  large  classes  of  terms  of  all  orders.     The  remain- 
ing terms  are  identical  with  corresponding  ones  in  the  expansion  for  the  true 
problem  (zero  couplings).     Although  still  of  all  orders,   they  have  a  sufficiently 
simple  structure  that  their  sum  represents  the  exact  solution  of  closed  integral 
equations. 

The  method  of  stochastic  models  is  introduced  in  the  present  paper  by 
application  to  a  linear  oscillator  whose  frequency  parameter  is  Gaussianly  dis- 
tributed over  Etfi  ensemble.     This  system  has  the  virtue  that  it  can  be  solved 
exactly.     Furthermore,   it  displays  great  sensitivity  to  inadequacies  in  approxi- 
mation schemes.     Neither  tiruncation  of  the  perturbation  series  nor  the  curaulant- 
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discard  approach  yields  admissible  approximations    (Section  2).     The  collective 
representation  and  the  general  model  are  formulated  in  Sections  3  and  h.     Expli- 
cit solutions  for  the  average  impulse-response  function  of  the  oscillator  then 
are  obtained  for  two  particular  models   (Sections  5  and  7).     They  suggest  the 
existence  of  a  sequence  of  model  solutions  which  converges  rapidly  to  the  axact 
solution  for  the  true  problem.     In  Section  8,  model   equations  are  obtained  for 
the  mean  and  covariance  of  the  anplitude  of  the  oscillator  when  driven  by  random 
forces.     The  generalization  to  non-Gaussian  frequency  distributions  is  described 
in  Section  9. 

In  Section  6,   approximations  for  the  average  response  function  are  examined 
which  represent  infinite  classes  of  terms  in  the  perturbation  expansion  for  the 
true  problem,  but  which  do  not  correspond  to  possible  stochastic  models.     Although 
they  are  very  plausible  in  terms  of  a  diagrammatic  representation  of  the  pertur- 
bation series,   these  approximations  have  pathological  characteristics.     This 
su^ests  that  correspondence  to  stochastic  models  may  be  a  useful  criterion  to 
help  judge  the  validity  of  partial  summations  of  perturbation  series  in  analogous 
situations. 

In  Sections  10  and  11,    stochastic  models  are  formulated  for  two  problems  of 
more  physical  interest:     the  Schrodinger  equation  of  a  particle  in  a  random 
potential  and  Burgers'  analog  to  turbulence  dynamics.     For  both  problems,   closed 
integral  equations  are  obtained  which  determine  the  average  Green's  function, 
the  amplitude  of  the  mean  field,   and  the  covariance  of  the  fluctuating  field. 
The  models  for  these  systems  have  an  intimate  formal  relation  to  those  for  the 
random  oscillator.     In  fact,   the  random  potential  problem  is   homologous  to  the 
oscillator  problem  in  the  sense  that  the  coupling  coefficients  characterizing 
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corresponding  models  are  identical  in  the  two  cases.  Many  results  for  the  random 
potential  problem  can  be  obtained  by  inspection  from  the  oscillator  results. 
A  comparison  of  the  model  equations  for  the  random  potential  and  turbulence  pro- 
blems illustrates  the  similarities  and  differences  involved  when  the  present 
method  is  applied  to  systans  which  are,  respectively,  linear  and  nonlinear  in 
the  dynamic  variables. 

In  a  paper  to  follow,  stochastic  models  are  formulated  for  classical  and 
quantized  nonlinear  oscillators.  Then  the  many-boson  problem  with  interparticle 
forces  is  treated.  This  problem  is  homologous  to  the  quantized  nonlinear  oscil- 
lator in  the  same  way  as  the  random  potential  problem  is  to  the  classical  random 
linear  oscillator.  Particular  attention  is  given  to  thermal  equilibrium.  The 
Einstein-Bose  distribution  law  is  derived  by  requiring  equilibrium  under  arbitrary 
infinitesimal  changes  in  the  coupling  among  systems  in  a  collection,  without 
assuming  the  grand  canonical  or  other  such  distribution. 

2.  THE  RANDOM  OSCILLATOR 

Let  the  amplitude  q(t)  of  a  linear  oscillator  satisfy 

dq(t)/dt  +  ibq(t)  -  0  ,  (2.1) 

where  b  is  a  real,  time-independent  parameter  which  is  statistically  distri- 
buted over  an  infinite  ensemble  of  realizations  of  the  oscillator.     We  shall 
be  interested  in  determining  the  function  G(t).  ■     (o  At)/  ,  where  {/  denotes 
ensemble  average  and  G     (t)  is  the  response  function  defined  for  "  oo  <  t  <  oo 

dGj.^(t)/dt  -  -  ibGj.^(t),  G^^(0)  -  1  .  (2.2) 
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00 

We  have,    imnediately,   G(t)  -     <(exp(-ibt)')    -  exp(-ibt)P(b)db,  where 

-00 

P(b)  is  the  normalized  probability  density  for  b«  Hence, 

0((o)  -  P(co)  ,  (2.3) 

where 

00 

G(co)  -  (2n)"  j  G(t)exp(iwt)dt  . 

-00 

Since  P(b)  >  0,  G(co)  must  satisfy  the  realizability  condition 

G(a))  -   |G(co)|.  (2.lj) 

A  particular  consequence  of  (2.U)  is 

|G(t)|  <  G(0)  -  1  ,  (2.5) 

which  also  follows  from  the  fact  that  Iq(t)|  is  a  constant  of  motion  in  each 
realization  of  the  oscillator. 

Now  suppose  that  P(b)  is  not  known  in  closed  fonn,  but  instead  is  speci- 
fied by  the  infinite  set  of  moments  /b),  <(b  ),  \h    ),  .  .  .  Then,  by  integrat- 
ing (2.2)  from  0  to  t,  iterating,  and  averaging,  we  may  generate  the  formal 
solution 

G(t)  -  1  ♦  X!  "^  (-i)''(b'')t"/ni  .  (2.6) 

n-1 

averaged 
Eq.    (2.6)  corresponds  precisely  to  the  perturbation  series  for  the^Green's 

function  in  certain  statistical  field  physics  problems.     Let  us  explore  its 

validity  for  the  present  problem  by  taking  the  example  of  Gaussian  P(b).     Then 


we  have 


G(co)   -  (2n  <b2  )  r^^^expi'  \  J/  fy^  })     ,  (2.7) 

G(t)  -  exp(-  |<b2)t2)  (2.8) 


But  let  us  suppose  that  we  do  not  know  this  closed  form  and  instead  are  given 
the  moment  values 

<b">-0       (nodd),  <b2")    -   (  |<b2))"2nj/nJ     .  (2.9) 

By  (2.6)  we  have 

G(t)  -  1  *  H  °°    (-  I  <b^  >t2)7nJ     ,  (2.10) 

n*! 

which,   of  course,  is  the  power  series  expansion  of   (2.8)   . 

The  following  observations  may  be  made  concerning   (2.10).     First,   it  is 
absolutely  convergent  for  all  t.     Second,  for  t  >  2/<b   >    '     the  convergence 
rapidly  becomes  very  poor  so  that  very  many  terms  must  be  taiken  to  obtain  a 
good  approximation.     Third,   if  the  series  is  truncated  after  any  finite  number 
of  terms,  we  have  G(t)  — >  oo,   t  — >  oo,  in  violation  of  the  basic  realizability 
condition  (2.5).     Thus,   at  no  finite  stage  of  the  iteration  treatment  do  we 
obtain  an  approximation  with  uniform  validity  for  all  t,  and,   in  pai'ticular, 
at  no  stage  does  the  spectral  density  G(co)  exist. 

Let  us  next  apply  a  second  approximation  scheme  which  has  been  widely 
used  in  statistical  field  physics.     From   (2.2)  we  may  obtain  the  infinite  set 
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of  coupled  equations 

dG(t)/dt  =  -i  <bG  (t))  G(0)  =-  1  ^ 

d  <bG^.^  (t)>  /dt  -  -i  <b^G^^(t)>         <bG^^(0)>  -  0  , 
d  (b^G  (t)>  M  -  -i  (b^G  (t)>        <b2G^^(0)>  =.<b^>. 

We  may  close  off  this  hierarchy  at  successively  higher  stages  by  taking  the 


(2.11) 


zeroth  approximation  that  b  and  G  (t)  are  statistically  independent,  and  then 
admitting  successively  higher-order  cumulants  of  the  joint  distribution 
(higher  'correlations'  in  the  language  of  statistical  field  physics).  Let  us 
again  assume  (2.9).  Then  the  appropriate  successive  closure  approximations 
are 

<bO|.^(t)>   -  <b>  G(t)  -  0, 


.2G^^(t)>  -  <b2) 


<b2G|.^(t)>  -  (b^)  G(t), 


VG^j(t)>  -  3<bO  <bG|.^(t)>^ 

(b'^G^^  (t)>  -  6  <b2>  (b^G  ^  (t)>  -  3  <b2>  2G(t),  (2.12) 

(Note  that  if  G  (t)  were  statistically  independent  of  b  then  all  these  rela- 
tions would  be  exact.)  Using  these  relations  in  turn  to  close  off  (2.11)  at 
successively  higher  stages,  we  obtain 

G(t)  -  1, 

G(t)  -  cos(b^t), 

G(t)    -  1+  5  C03(/Tb^t), 

1  r  1/2     1      1  r  1/2     -I 

G(t)  -  3(3+/5)co3l  (3-/5)      b^t     +  ^(3-/5)008 1  (3 +y5)      b^tl  (2.13) 
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Where  b^  -  (b^)  ^/^ 

Beyond  the  zeroth  stage,  which  yields  identical  results  in  the  two  cases, 
the  sequence  (2.13)  is  distinctly  superior  to  the  approximations  obtained  by 
truncating  (2.10).  All  members  of  the  sequence  satisfy  (2.'i)  and  (2.5).  In 
cormon  with  the  iteration  scheme,  the  first  n  even  derivatives  of  G(t)  at  t  =  0 
are  correct  in  the  nth  approximation.  However,  there  still  is  no  uniform 
validity  in  the  sense  G(t)  — >  0,  t  — >   oo  .  None  of  the  moments 

00 

]     t"G(t)dt     (n  »  1,2,...)  exist  for  any  approximation  in  the  sequence,  whereas 
-'O 

they  all  do  for  the  exact  solution.     Alternatively,   we  may  note  that  G(i:o),  which 
is  smooth  in  the  exact   solution,    is  a  sun  of  6  functions  in  any  of  the  cumulant- 
discard  approximations.     The  convergence  to  the  exact  G(t)  is  still  very  poor 
for  t  >  2A^  . 

It  is  clear  that  the  random  oscillator  exhibits  in  acute  form  certain 
shortcomings  of  the  iteration   (perturbation)  and  cumulant-discard  approaches 
to  dynarriical  equations  which  are  nonlinear  in  stochastic  quantities.     Both  for 
this  reason  and  because  of  its  simplicity,  we  shall  use  the  random  oscillator  to 
illustrate  the  alternative  approach  which  is  the  subject  of  this  paper.     The 
sensitivity  to  inadequacies  in  the  method  of  approximation  arises  because  the 
solution  to  the   'unperturbed'   equation,   obtained  by  replacing  the  right  side  of 

(2.2)  with  zero,   has  a  monochromatic  spectrum.     In  this  respect  it  resembles 
certain  limiting  cases  of  statistical  field  theory  problems  which  are  of  current 
interest  and  to  which  our  approach  will  be  applicable.     Examples  are 

a       quantum-nechanical  particle  in  a  random  potential  in  the  WKBJ  limit,   tur- 
bulence at  infinite  Reynolds  number,   and,    in  a  less  direct  sense,   a  second- 
quantized  many-boson  system  at  very  low  tanperature. 
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3.      COLLECTIVE  REPRESENTATION  FOR  A  SET  OF  OSCILLATORS 

We  shall  now  describe     a        dynamical  representation  which  is  appropri- 
ate for  formulating  the  stochastic  models  promised  in  Section  1,     We  shall 
introduce  the  representation  formally  and  then  give  a  physical  interpretation 
and  a  comparison  with  more  familiar  concepts. 

In  Section  2  we  treated  an  ensemble  of  realizations  of  a  single  oscillator. 
Now  let  us  consider  a  collection  of  M  oscillators   (M  =  23  +  1,   S  =  positive 
integer)  whose  frequencies  are  identically  and  independently  distributed  over 
an  ensemble  of  realizations  of  the  collection.     We  shall  be  interested  in  the 
limit  M  — >  00,    so  that  in  reality  we  are  introducing  a  two-dimensional  dis- 
tribution,  so  to  speak.     In  place  of  the  frequencies  and  amplitudes  of  the 
M  individual  oscillators,   let  us  adopt  the  collective  parameters  and  coordi- 
nates 

b„  -  M'-^/^  y    exp(i2nan/M)br-i  , 
q^(t)  -  M'^/^  Y1     exp(i2nan/M)qj.     (t),  (3.1) 

(a  -  0,  +  1,    ...,    +  S)  (n  -         1,  2,    ...,  M)     , 

where  b  ^ -,   and  q  ;„■](*)  are  the  frequency  and  amplitude  of  the  nth  oscillator. 
The  identities 

M"^  ^    expri27ia(n-m)AJ    "  6^  „  ,       ^'^  H    exp|i2nn(a-p)/M]   -  5^  (3.2) 
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yield 

-  M-1/2 


bn   -  VT^'^  y    exp(-i2TOnA)b    , 
qr-,(t)  -  M"^/^  ^     exp(-i2nan/^)q^(t)    .  (3.3) 


Let  us  adopt  hereafter  the  cyclic  convention  o  +  M  ;  a  ,   which  clearly  is  con- 
sistent with  (3.1)  and  (3.3)   . 

From   (3.1)  and   (3.2)  we  easily  find 

b       ■  b  *  7~  b  b       .  F  br  -,^  , 

-a         a  '  Z_^  ,  -a  -  ^^  [n]      » 

where  we  have  used  bi-  -i  -  bp-,   .  From  the  equations  of  motion 

dqr  -,(t)/dt  ♦  ibr„-iqrn-j(t)  -  0  ,  (3.5) 


we  find 


dq^(t)/dt  -  -iM'^/2  ^  ^bVs^*^  *  ^^*^^ 


where  a  -  p  is  to  be  interpreted  according  to  the  cyclic  convention.  This  shows 
that  the  new  coordinates,  in  contrast  to  the  old,  are  dynamically  coupled.  By 
(3.'4),  YL    ^J^'^'^     (*)  ^^  ^  constant  of  motion.  Let  Op  -i  (t)  denote  the  solution 

of  (3.5)  with  qn(0)  -  5    and  let  Q   (t)  denote  the  solution  of  (3.6)  with 
[rj      r,m        a,Y 

q  (0)  ■  5    .  Since  the  individual  oscillators  are  uncoupled,  we  have 
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0^      ^(t)  -  5       G-     -,(t).     Hence,  by  (3.1),    (3.3),   and  the  linearity  of  the 
[^njHij  n^ni    in^nj 

equations  of  motion,  the  response  matrices  in  the  two  representations  are 
related  by 

O^^^it)  '  }f^J2    exp[i2n(a-Y)n/^]G  (t).  (3.7) 

The  statistical  properties  of  the  b  are  easily  found  from  the  assumption 
that  the  b  j-  -.  are  identically  and  independently  distributed.  By  (3.1)  and 
(3.2),  we  have  immediately 

(b^bp  .  .  .  b^N^  -  0      (o  +  p  +  .  .  .  ♦  ti  /  0).      (3.8) 

For  Gaussian  b_  .,   all  odd-order  moments  vanish,   and  we  have 

(br  .bp  -,>     -  5         (b^>, 
[n]    [m]  /  n,m  ^     ^  ' 

<bnbnbr.bn>       -(5       6         +6       6         +6       5       )(b^>^,  (3.9) 

[nj    [mj    [r]    [s]/  n,m  r,8         n,r  m,3         n,s  m,r     ^  '  \^'^/ 

where     (b  )     -  /bj-  -,    )  ,  which  is  the  same  for  all  n.     Then,  by   (3.1), 

/b     ...b\-0  (odd  number  of  factors). 


<  Vp>     ■  ^a+p  ^^^  (^a+p  -  -0,a+p 


2\2 


/b  b„b  b    N    ■   (5     „5         +  5       5„       +  6       6„      )  (b    )  (3.10) 

\    o  p  Y  e'^         ^  o+p  Y+e         a+Y  P+e         o+e  p+Y  ^     >  \J'^^J 

For  any  univariate  bj  -,   distribution  we  find  from  (3.7)  . 


-  lU  - 

<O^^Y(t)G*^^(t)>  -  5^^^|G(t)|2  *  M-^[<|G^j(t)|2>-iG(t)|2]  •         (3.11) 

Here  G  (t)  denotes  G^   -,  (t),  which  is  the  same  for  all  n,  and  G(t)  -  (G  (t))  , 
as  in  Section  2.  In  the  limit  M  — >   co,  with  which  we  are  concerned,  (3.11) 
gives 

<Q^^^(t)>  -  G(t),        <|G^^^(t)  -  G(t)|2>  -0(M-^), 

(0^^^(t)G*^^^(t))  -0(M-^),    <Va-Y,a^*^>  -  0(M-^/2j   (a  /  y)  •  (3.12) 

These  relations  show  that  the  variance  of  G       (t)  «  M~     Y~    (L     _  (t)  vanishes 

a,  a  ^  ]n,nj 

in  the  limit.     That  is,   G       (t)  is  statistically  sharp.     They  fvirther  imply  that 
the  effective  dynamical  coupling  between  any  given  pair  of  degrees  of  freedom 
q     and  q     is  infinitely  weak  in  the  limit.     Eqs.    (3.12)  were  obtained  without 
explicit  reference  to  (3.6)  but  their  (^namical  implications  may  also  be  inferred 
from  the  latter.     The  direct  dynamical  coupling  of  q     to  any  o     arises  from 

only  one  of  the  M  terms,   each  0(M"  '    ),   on  the  right  side  of   (3.6).     Eqs.    (3.12) 

-1/2 
show  that  the  effective  coupling  is  still  0(M     '    )  when  the  indirect  interaction 

of  q    and  q.^  through  all  the  other  degrees  of.  freedom  is  included.     The  absence 

of  fluctuations  in  G       (t)  in  the  limit  is  consistent  with  the  fact  that  this 
a,a 

f\mction  is  determined  by  the  simultaneous  interaction  of  q^  with  an  infinite 
number  of  other  degrees  of  freedom;  q     exhibits  negligible  self-coupling,  in 
contrast  to  q^  ,    • 
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It  is  apparent,  both  from   (3.1)  and  the  convolution  structure  of   (3.6), 

that  the  q     have  a  close  formal  relation  to  Fourier  coefficients.     The 
a 

physical  significance  of  the  new  representation  is  best  brought  out,   in  fact, 
by  a  comparison  with  analysis  into  wavenumber  or  frequency  components.     Let 
)^(x,t)  be  a  scalar  field,   associated  with  an  extended  (one-dimensional)  dynami- 
cal system,  which  is  described  by  an  ensemble  statistically  invariant  under 
translation  (e.g.,  (}^(x,t)/'(x',t))    •  ( ;^(x+y,t))^(x'+y,t )>     for  all  y)   .     The 
natural  coordinates  for  describing  the  field  then  are  wavenunfcer  components, 
which  change  only  by  a  phase  factor  under  translation.     Suppose,  instead,    the 
ensemble  were  invariant  under  time  displacement.     Then  the  natural  coordinates 
would  be  frequency  components. 

Physical  systems  usually  are  neither  statistically  homogeneous  nor  station- 
ary.    However,  if  we  form  a  collection  of  identically  distributed  individual 
systems,  then  obviously  (and,  it  will  appear  at  first  sight,  trivially)  there 
is  statistical  invariance  under  permutation  within  the  collection.     The  new 
representation  is  natural  in  the  presence  of  this  invariance  in  the  same  way 
that  a  wavenui*er  representation  is  natural  when  there  is  translational  invari- 
ance* 


To  examine  the  analogy  further,  let  us  take 

L/2 

y^j^(t)  -  L"^   f       5^(x,t)exp(-ikx)dx  (k-2naA-,  a-0,+1,...), 

U/2 

where  we  adopt  the  customary  device  of  making  the  field  cyclic  with  a  period  L 
which  is  as  large  as  we  wish  compared  to  any  relevant  correlation  length.  Let 
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us  divide  L  into  very  many  segments,   each  still  very  large  compared  to  any 
correlation  length.     Then  each  segment  contains  a  sub-system  which  has  only 
a  negligible  statistical  dependence  on  its  neighbors.     Furthermore,   it  is 
plausible  to  suppose  that   (over  times  which  are  not  too  large)  each  sub-system 
has  only  a  negligible  dynamical  interaction  with  its  neighbors.     Then  we  validly 
may  regard  the  set  of  sub-systems  as  analogous  to  the  collection  of  perfectly 
independent  systems  used  above  in  defining  the  q  .     Considered  in  this  way, 
the^'.    and  the  q     (for  large  M)  play  essentially  similar  roles.     Both  are  linear 
combinations  of  the  physical  coordinates  of  a  very  large  number  of  effectively 
independent  systems. 


U.     FORMULATION  OF  MODEL  PROBLEMS 

Consider,  instead  of  (3.6),  the  more  general  equations 


dqjt)/dt  -  -iM-1/2  r    5fa,p,a.^^pVp(*^     '  ^^'^^ 


whei^  ^    a       Q  is  independent  of  t  and  the  same  for  every  realization  in  the 
ensemble.  We  shall  be  interested  in  stochastic  assignments  of  0„  „  „  _  in  the 

sense  that  this  quantity  will  exhibit  random  changes  in  value  as  a  and  p  are 

changed.  By  (3.3)  and  (3.2),  (U.l)  implies 


Where 


Thus,  the  individual  oscillators  in  the  collection  now  are  dynamically  coupled. 
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When  ^  „   .  -  1  for  all  a  and  S,  then  Ap    ->  ■  5   6    so  that  we  recover 
a,p,a-p  ^'       |_n,r,  sj    r,  n  s,n 

the  original  collection  of  uncoupled  oscillators.  The  quantities  q^ -^  (t  )qp  ..  (t) 
are  no  longer  constants  of  motion  in  the  general  case.  However,  we  shall 
require 


Then,  since  b  „  ■  b.  ,  we  find 
-p    P 

n  '--''-  -■  a 

The  response  matrix  corresponding  to  (U.l)  satisfies 

dG   (t)/dt  -  -iM-^/2  J-     ^  b  G     (t)  ,     G   (0)  -  5    .   (U.6) 

a,Y  ' — ^g  °-*Pf^'P   P  <i"p*Y  ijY      a>Y 

Suppose  that  we  carj*y  out  a  formal  iteration  expansion  of  (U.6).  The  coefficient 
of  t  in  the  resulting  power  series  for  G   (t)  is  a  sum  over  products  of  n 
factors  ^  and  n  factors  b.  It  is  clear,  from  the  initial  condition  and  the 
w^  the  Indices  combine,  that  in  each  product  the  sum  of  the  indices  of  the  b 
factors  must  be  a-y  •  Hence,  by  (3.6), 

(G^^^Ct))  -0        (a  fir)  (U.7) 

By  (3.1),  (3.3),  and  (U.?),  we  have 

<G     (t)>  -  M~^  ^  exp[-i2n(n-m)aAl]  (G^^^(t))  .        (^.8) 


18 


Therefore,  if  (g   (t))  is  independent  of  a,  we  have,  by  (3*2), 

G.      it)  -  5^     G(t)  ,  0(t)  -    <0       (t)>     ,  (U.9) 

Ln,ifij  n,it>  ^   a,  a       ' 

where  G(t)  has  now  the  same  meaning  as  in  Section  3«  We  shall  be  concerned 
only  with  0  assignments  which  yield  (U.9),  and  presently  we  sYsb.11   exhibit  cer- 
tain of  their  properties. 

We  wish  now  to  develop  an  expression  for  d  (G   (t))  /dt  by  expanding  the 
right  side  of  (h.6),  with  y  ="  Q  .  Writing 

H  .  p  (t)  -iM-^/2  V  ,^j  (^^3Q) 

a-p,p,a  a,p,a-p  (3  a-p,a 

we  may  develop     (H     .  „      (t))     in  powers  of  t  by  using  the  iteration  expansion 
ci-p,p,a 

r-l 
for  G     „     (t).     The  coefficient  of  t         in  the  resulting  series  is  a  sum,   over 

^-P**'   all 

the  indices  of  b  factors  except  b^,  of  products  of  r  factors  b  and  r  factors  0 
A  p 

fincluding  the  factors  b-  and  0  „   „  which  appear  explicitly  in  (U.IO)] .  Let 
<-  p  a,p,a-p 

us  take  Gaussian  b,  ^,  (We  shall  return  to  the  general  case  in  Section  9.) 

CL^td   (3.  ui;L"J 

Then,  by  (3.8),  only  the  odd  powers  of  t  survive  and  in  the  coefficients  of 
these  the  indices  of  the  b  factors  must  be  equal  and  opposite  in  pairs.  In 
this  way  we  find 


<Vp,p,e^*^>  '^'^iT^    L  (-l)\„.p(a,?,a-p)<b2>V"-V(2n-l)i,   (U.ll) 


where  m"  C^  ^  (a,p,a-p)  is  a  sum  of  products  of  2n  factors  ^  and  the  index  p 
2njp 

(p»l,2,...,2nV2"nJ)  labels  (in  an  arbitrary  order)  the  contributions  which 
arise  from  all  the  possible  pairings  of  the  b  factors. 
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The  Cp  .  (a,p,a-p)  through  n  «  2  are 

%}1^    )  "  «'^  E^  ^a,p,a.p  Cp,-p,a  S^a,Y,a.Y  ^a-Y,-Y,a  , 

%.2^'>     '     ""^II^52f^,p,a-P^a-p,Y,a-p.Y^a-p-Y,-Y,a-p^a-p,-p,a, 

Cl,.3(  )  -  ""^L^^a,p,a-p  Cp,Y,a-p-Y^a.p-Y,-p,a-YCY,-Y,a       (U.12) 

The  (^  factors  in  (U.12)  other  than  ^  _   _  are  written  from  the  right  in  the 
order  they  arise  in  the  iteration  process.   (The  values  of  p  are  assigned  arbi- 
trarily. ) 

Each  C_   (a,p,a-p)  is  closed  in  the  sense  that  the  initial  index  on  any 
<injp 

factor  equals  the  final  index  on  the  factor  to  its  immediate  left  [when  ordered 

as  in  (b.l2)]  while  the  middle  indices  are  equal  and  opposite  in  pairs.  This 

permits  a  systematic  diagrammatic  representation  of  the  C   „(a,P,a-p)  and 

<:n,  p 

therefore  of  (U.]!).  With  each  ^  ,   or  5*  .   let  us  associate  a  vertex  as 
shown  in  Figure  1.  The  Cp  .  (a,p,a-p)  then  may  be  obtained  by  the  following 
rules: 

Connect  2n  points  by  2n  solid  line  segments  to  form  a  single  closed  loop; 
then  connect  all  the  points  in  pairs  by  n  dashed  line  segments  to  form 
a  closed  diagram  of  2n  vertices.  Equip  all  the  solid  line  segments 
with  arrows  pointing  in  the  same  sense.  Choose  one  vertex,  equip  its 
dashed  line  segment  with  an  ingoing  arrow,  and  label  its  three  line  seg- 
ments to  correspond  to  0  „   o  Iri  ^he  sense  of  Fig.  la.  Call  this  the 
fixed  vertex  and  identify  it  by  circling.  Label  the  remaining  dashed 
line  segments  Y>  «>  ...  in  any  order  and  equip  them  with  arrows  (>rtiose 
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direction  does  not  matter     ).     Complete  the  labelling  of  the  solid  line 
segments  so  that  the  sum  of  the  indices  labelling  the  ingoing  lines 
equals  the  sum  of  those  labelling  the  outgoing  lines  at  every  vertex. 
Now  write  the  product  of  the  j2f  factors  associated  with  all  the  vertices 
according  to  Fig.  la,   or  according  to  Fig.  lb  and   (U.U).     For  each  n 
there  are  2ni/2"nJ  distinct  diagrams  of  this  type,   corresponding  to 
the  2nJ/2"ni  ways  of  connecting  the  vertices  by  dashed  lines  after  one 
is  chosen  as  the  fixed  vertex.     Each  diagram  corresponds  to  one 
C     ^    (a,p,o-p)    (according  to  an  arbitrary  rule  for  assigning  the  values 
of  p)  and  the  latter  is  equal  to  ir"^  times  the  sum  over  r,t,    ...  of 
the  associated  ^  product. 

The  diagrams  associated  with   (U.12)  are  shown  in  Fig.   2. 

It  is  clear  from  (U.ll)  that    (G    jj(t))     will  be  independent  of  a,   and 

therefore   (h.9)  will  hold,   if 


M-'r     %p(a,p,a-p).C2„.p,  (U.13) 


where  Cp    is  independent  of  a.  In  this  case, 

G(t)  -  1  +  e"*  E  (-l)\n-o  <^^  >"t^"/2n'  •  (i»-l»*) 

n«l   p        *^ 

We  shall  be  concerned  hereafter  only  with  ^  assignments  such  that  (U.13)  is 
satisfied  when  M  — >  oo  .  The' C.    may  be  interpreted  as  moments  of  the  distri- 
bution of  the  quantity  ^  over  the  set  of  index  values  u.  and  X. 
^i,  A.,  |A-  A. 

Let  us  associate  with  Cg^   the  diagram  for  C^^^^{a,^,a'^),   but  with 
index  labels  and  dashed  line  arrows  (which  are  now  superfluous)  omitted.  By 
(U.13)  we  have  C^^^^  -  M*^  ^    Snip^^''^**'"^^  '  Recalling  the  cyclic  coiiven- 
tion  11  -  ii  +  M  (any  ii)»  we  see  that  the  sunmation  in  this  expression  is  equiva- 
lent to  one  over  all  M  values  of  all  the  indices  labelling  lines  in  the  diagram, 
subject  only  to  the  sum  coiviition  at  each  vertex.  Consequently,  the  expression 
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is  independent  of  which  is  the  fixed  vertex;  its  value  depends  only  on  the 
order  and  topology  of  the  diagram. 

It  is  important  to  point  out  that  (3.12)  is  valid  for  general  (^'s  satis- 
fying (U.13).  In  particular,  G   (t)  is  statistically  sharp  (M  — >  oo).  In 
the  original  case  (all  0's  unity)  this  was  so  because  q  interacted  simul- 
taneously with  all  the  other  degrees  of  freedom,  and  negligibly  with  itself. 
These  properties  clearly  characterize  the  general  case  also,  provided  the  0's 
are  bounded  as  M  — >   oo  .  The  validity  of  (3.12)  in  the  general  case  is  easily 
demonstrated   for  any  power  of  t  in  the  iteration  expansions  of  the  left  sides 
of  the  equations,  if  one  uses  (3.10)  sind  the  fact  that  the  0's  are  the  same  for 

all  realizations.  We  shall  use  the  abbreviations  C  and  C(  )  to  denote  C^    and 

2njp 

C„   (a,p,a-p),  respectively,  when  it  is  not  desired  to  specify  particular  sub- 
^njp 

scripts  and  arguments. 

Let  us  define  a  reducible  C  as  one  which  may  be  factored  into  two  or  more 
C's  of  lower  order,  and  an  irreducible  C  as  one  which  may  not.  Let  us  define 
a  reducible  C(  )  as  one  which  may  be  factored  into  the  product  of  a  lower-order 
C(  )  with  one  or  more  C's,  and  an  irreducible  C(  )  as  one  which  may  not.  It 
follows  that  each  reducible  C  is  a  product  of  irreducible  C's  and  each  reducible 
C(  )  is  a  product  of  an  irreducible  C(  )  with  irreducible  C's.  It  is  easy  to 
see  from  our  rules  that  reducible  C(  )'s  and  C's  (and  only  they)  are  associated 
with  diagrams  in  which  there  is  a  part,  or  parts,  connected  to  the  rest  of  the 
diagram  by  only  solid  lines.  Thus,  C^  ,  «nd  C,  _  are  irreducible,  but  C,  ,  and 
Cj^,2  are  reducible.  Using  (U.13)  we  find  C^  .  -  Cj^  _  ■  (^2il^^  * 

Let  us  write  each  C(  )  which  appears  in  (U.ll)  as  the  product  of  an  in*e- 
ducible  C(  )  and  irreducible  C's,  and  then  collect  all  the  terms  proportional 
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to  each  irreducible  C(  ).  We  obtain  a  result  of  the  form 

< Vf,p,a(*'>  ■  ""'  E^  E  '"   (-1)"  %p(°.P.-P"^n,p('>  '      *•«' 

irr 
where  ^    denotea  the  stm  over  irreducible  diagrams  only.  The  ^„  .r,^*^ 

depend  on  the  values  of  the  irreducible  C's  but  are  independent  of  a  and  p. 

Each  ^  .  (t)  contains  all  (odd)  powers  of  t  which  are  >  2n-l,  since  each 

C    (a,p,a-p)  appears  in  reducible  C(  )  »s  of  all  orders  >  2n.  The  J^^  „(*) 

turn  out  to  have  simple  expressions  in  terms  of  ^b  /  and  G(t)  which  may  be 

found  by  comparing  the  explicit  power  series  for  ^  (t)  and  G(t).  However, 

the  same  result  may  be  obtained  more  transparently  by  a  variational  procedure 

which  provides  certain  dynamical  insights. 

There  are  M(M-l)  0's,  and  only  M  sums  Yi    '^2n-p^°*^*°~^^  ^^^   given  n  and  p. 

In  the  limit  M  — >  OD  it  will  be  possible,  therefore,  to  make  wide  classes  of 

variations  A  ^  such  that  (U.h)  and  (U.13)  continue  to  hold  and  such  that 

A  C„_   -  0  for  all  finite  n.  Under  these  constraints,  /\i^     ^(t)  ■  0.  Con- 
2n;p  <injp 

sequent ly,  we  have 


■I       irr 


Now  consider  a   (finite)  variation  A  gf    .     «  „  -  A(2^,   „  „_q  for  a  particular 

a  and  p,  with  all  the  other  ^f's  fixed.     We  may  vary  the  real  and  imaginary 

parts  of  of    -     -       independently.     Identical  results  are  obtained  by  supposing 
a-p,-p,a 

^„  D     o  «  *o  ^^  whil*  ^«  ft     ft  n  "  ^n  ft  n  ft  is  held  fixed,  and  we  shall 
a-p,-p,o  a-p,-p,a         a,p,a-p 

adopt  the  latter  procedure.     Then,  by  (U.12)  , 

AC2,i(a,p,a.p)  -  Sf^^^^^.p  ^Cp,.p,a  ,  (U.17) 
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while  from  (b.l3)  we  see  that  A  C-    -  0(M"  )  for  all  n  and  p.  Thus  the  con- 

straints  stated  above  are  satisfied  for  M  — >  cd  .  By  (Li.6),  G  .   (t) 

a-p,a 

satisfies 

The  effect  of  the  variation  A^    r,     a       is  to  produce  on  the  right   side  of 
(U«18)  the  additional  terra 

-iM'^/^  Acf    „     „     b  „G       (t) 
a-p,-p,a  -p  a,a^ 

-1/2 
which,  we  note,  is  0(M  '  ).  Now  we  recall  that  G  „  (t)  is  simply  the  ampli- 


tude q  _o(t)  \mder  a  particular  initial  condition  at  t  ■  0.  Therefore,  to 
order  M  '  ,  we  have 

t 
^°a-p,a(*)  -  j  Vp,a-p(*-^)h""^^'^Cp,-p,aVa,a(^^]^  '         ^^-^^^ 

since  G     _       o  is  "the  diagonal  response  function  for  q     _  and  the  perturbation 
a— p,a-p  a— p 

does  not  affect  the  initial  condition.      [Note  that  AG       (t)  and  AG     »       o(t) 

o,fli  o—p, u—p 

are  0(M~  )  under  our  constraints.]     Referring  to   (U.IO),  we  obtain 


^     a  a     (t))     immediately.     It  is  clear  that  our  variation  gives 
ACg       (a,p,o-p)  -  0(M     )  for  all  irreducible  diagrams  with  n  >  1.     Then  by 
(U.16)  and  (U.17)  we  find 

t 

^;l^*>  "  [    <V-pVp,a-3(*-'^^a,a(«>>  ^  ('^•20) 

in  the  limit  M  — >  00.  As  we  have  noted  previously,  0   (t)  and  0  .   „(t) 
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are  statistically  sharp  in  the  limit.  Therefore,  by  (U.9)  and  (3.10), 

t 
^.^(t)  -  (b^)  f  G(t-s)G(s)ds  .  (U.21) 

The  higher  C>  .  may  be  found  by  similar  analysis  based  on  more  general 
variations.  The  result  is 

^n,p-  <b^>''G('I^G)^"'^  »  (i^-22) 

where  G(>(cG)     is  a  repeated  convolution;  e.g.  (for  argument  t), 

t     s     s' 
G(>j^G)^  -  [  ds  f  ds'  f   ds"  G(t-s)G(s-a  •)G(s'-s"  )G(s"  )  . 
^0     0     0 


Collecting  the  appropriate  relations,  we  have  the  final  result 
00         irr  ^  ^ ,  „  «„  1 

dG/dt  -  r    r    (-i)\„.n  <^  >  <5(  *Q)     ,      G(o)  - 1 

n-1       p  ^""'P 

The  value  of  this  infinite-series,  integro-differential  equation  for  G(t)  is  that 


(U.23) 


only  the  irreducible  C-    appear  explicitly. 


5,  THE  RANDOM  COUPLING  MODEL 

We  shall  now  consider  a  particular  stochastic  assignment  of  the  $('s.  Let 

Where  9  «   „  is  real  and  satisfies 

e  „  „    -  -©       .  (5.2) 

a-p,-p,a      a,p,a-p 
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For  each  choice  of  a  and  S,  let  ©or.  take  a  value  at  random  in  the  interval 

0  to  2n,  subject  only  to  (5.2).  The  value  must  be  the  same,  of  course,  for 

every  realization  in  the  ensanble.  Now  let  M  — 1>   oo.  Clearly  this  assignment 

satisfies  (U.U).  In  addition,  it  jrields  10  -   o  I  "  1»  and  therefore  retains 

unaltered  the  strengths  of  the  individual  dynamical  couplings  of  pairs  q  ,  q^_- 

which  characterize  (3.6).  Now,  however,  the  phages  of  the  couplings  are  completely 

unrelated  for  different  pairs.  We  shall  call  the  present  choice  the  random 

coupling  model. 

Referring  to  (U.12)  and  (U.13),  we  find  C^.^  -  1,  as  in  the  true  problem 

(all  0's  unity).  Consider  C.  .,  however.  Each  product  in  the  sum  has  modulus  1, 

but  the  phase  of  the  product  changes  at  random  with  p  and  y.     Consequently, 

C,  -  ■  0(M"  )♦  In  a  similar  fashion  we  see  that  the  only  C„    which  survive 
U»3  2njp 

in  the  limit  are  those  in  which  the  product  of  ff's  consists  entirely  of  conjugate 
pairs  and  which  therefore  are  reducible  to  powers  of  C_  ^ .  Consequently,  for 
M  — >  00  . 


C-    -  0    (n  >  1)   (all  irreducible  diagrams).       (5.3) 
<in;p 

It  follows  that  (U.23)  reduces  to  the  closed  form 

dG/dt  +  <b^>  G*G  -  0  ,      G(0)  -  1  .  (5.1l) 

Eq.    (5.U)  is  readily  solved  by  Laplace  transformation.     We  find 

G(t)  -  J^(2b^t)A^t  ,  (5.5) 

fn-^[l  -   (co/2b^)2]l/2  (|co|<  2b^,) 

^^"^■|0  (|co!>2b^),  (5.6) 


-  26  - 

where  b^  -  (b  )  '     as  before.  Considered  as  an  approximation  to  (2.8)  and 
(2,7),  the  present  results  display  a  type  of  uniform  validity  which  is  absent 
in  any  finite  stage  of  the  iteration  or  cumulant-discard  schemes  discussed  in 
Section  2.  The  spectral  density  (5.6)  is  continuous,  and  (2.U)  and  (2.5)  are 
satisfied.  All  the  derivatives  of  G(a>)  exist  at  co  ■  0  so  that  all  the  moments 

f  t'^G(t)dt  exist.  The  results  (2.8)  and  (5.5)  are  compared  in  Fig.  3. 

0 

It  is  important  that  certain  of  these  properties  could  have  been  predicted 

from  the  sole  fact  that  (5.!i)  is  an  exact  equation  (M  — >  oo)  for  a  realizable 

0  assignment  satisfying  (U.U)  and,  hence,  for  a  conservative  dynamical  problem. 

We  recall  that  G^  ^(t)  -  q^(t)  when  q^(0)  -  5^  .  But  then  ^  q^(0)q^*(0)  -  1, 

and,  since  ^I'lv'^*^^  ^*^  ^^  ^  constant  of  motion,  we  have  |G   (t)|  <  1,  whence 

(2,5)  readily  follows.  To  establish  (2.I4),  let  us  make,  for  each  realization, 

a  similarity  transformation 

(B   independent  of  t)  such  that  (h.l)  is  brought  to  the  diagonal  form 

dq^'(t)/dt  +  la3^q^»(t)  -  0  .  (5.7) 

Since  r~  q^'q.^'     "  ^  q  q       is  a  constant  of  motion,  the  co     are  real.     Now 

if  G       »(t)  is  the  response  matrix  of  the  new  variables,   G^  ,'(t)  •  5^  ^exp(-ico  t), 
Y»  8  T>  6  T,  6  Y 

and  therefore 

Hence  0       (co)  is  real  and  non-negative  in  each  realization,  which  implies  (2.U). 
Finally,  we  note  that  the  model  problem  resulting  from  a  general,   realizable  j( 
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assignment  involves  the  interaction  of  an  infinite  number  of  degrees  of  freedom 

q  when  M  — >  oo .  From  this,  and  the  fact  that  the  b  are  continuously  dis- 
a  CL 

tributed  over  the  ensemble,  we  may  anticipate  that  G(co)  exhibits  a  continuous 

or  band  structure  and  that  G(t)  — >  0,  t  — >   oo  . 

It  is  possible  to  understand  qualitatively  why  the  complex  detailed  dynamics 

of  the  random  coupling  model  lead  to  a  simplification  of  the  statistical  dynamics 

and  to  closure  of  (U.23).  The  function  (G   (t))  describes  the  decay  of  q^ 

^  a,G   '  0. 

due  to  transfer  of  an  initial  excitation  q  (0)  -  1  to  the  rest  of  the  degrees 

a 

of  freedom.  In  general,  the  decay   requires  that  certain  phase  relations  be 

set  up  between  q  and  the  other  amplitudes.  A  phase  relation  between  q  and 

q  .  can  arise  either  from  direct  dynamical  coupling  (involving  the  coupling 
a-p 

coefficients  -iK     '     0    „   „b_  and  -iM~  '0  „  „  b  .)  or  from  indirect  coup- 
a,p,o-p  p  a-p,-p,a  -p 

ling  through  chains  of  other  modes  q^.  In  fact,  each  term  in  the  irreducible 

diagram  expansion  (U.l^)  may  be  regarded  as  describing  the  transfer  of  excitation 

from  q  to  q  -  along  the  chain  of  intermediate  modes  represented  by  the  directed 
a.  a"p 

solid  line  segments  in  the  associated  diagram.  The  closing  of  the  solid  line  on 

itself  then  represents  the  reaction  on  mode  a  and  the  consequent  diminuation  of 

q  .  The  factors  G,  >*iose  repeated  convolution  yields  the  ^  .  ,  incorporate  the 

effect  of  the  dynamical  interaction  as  a  whole  on  the  transfer  process.  This 

effect  is  to  relax  the  phase  relations  set  up  along  the  chain. 

In  the  random  coupling  model,  only  the  direct  interaction,  associated  with 

C_  -  ,  is  effective  in  the  transfer  of  excitation.  The  contributions  associated 

with  the  indirect  paths  of  interaction  cancel,  when  summed  over  all  possible 

intermediate  modes,  because  of  the  randan  phases  of  the  0's.  The  coupling  of 

q  and  q  „  to  all  the  rest  of  the  modes  therefore  affects  ( H  .  -  (t)  )    only 
a,  a— p  ^  i-p,p,a   ' 

by  relaxing  the  phase  relations  induced  by  the  direct  interaction  of  these  two 
modes. 
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All  C^        which  are  expressible  as  powers  of  C„  ^   have  the  value  unity  in 
2njp  d;i. 

the  randOTi  coupling  model,   and  all  other  C^„.„  vanish.     Thus  we  see  from   (U.lU) 

that  the  power  series  for  the  model  G(t)  consists  of  a  particular  sub-set  of 

terms  of  all  orders  from  the  corresponding  series  for  the  true  problem  (all 

C-         unity).     The  terms  retained  are  all  those  whose  associated  diagrams  can 
Znjp 

be  reduced  to  that  for  Cp  ^    (Fig.   2)  by  iterating  any  number  of  times,   on  any 
solid  lines,  the  contraction  operation  shown  in  Fig.   Ua.     Examples  of  included 
diagrams  are  shown  in  Fig.   Hb.      It  follows  readily  from  (5.5)  that  the  number  of 
diagrams  of  this  type  with  2n  vertices  is 


(-1)' 


^2n 


(dt)' 


Jl(2t) 


t  -  0 


6.     INAEMISSIBLE  HIGHER  APPRCfXIMATIONS 

For  the  true  problem  (all  Cg^      -  l),    (U.23)  reads 


dO/dt  -  21°°    (-!)%„  <b2>  '^GOfG)^"-^  ,  0(0)  -  1  ,  (6.1) 


n=l 


where  S-  is  the  number  of  irreducible  diagrams  of  order  2n.  The  first  few 
2n 

S-  are  S„  -  1,  S,  -  1,  S,  «  U,  Sg  -  2?.    The  relative  success  of  (5.1«)  as 

an  approximation  to  (6,1)  suggests  that  we  seek  higher  approximations  satisfying 

equations  of  the  form 


vn„   /,.2\  n„,,,.  „v2n-l 


dG/dt-^   (-1)  S   <bO  G(*0)'="-^  ,     Q(0)  -  1   (R>1),      (6.2) 
n-1 


which  we  obtain  by  giving  all  irreducible  Cg^^   the  value  one,  n  <  R,  and  the 
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value  zero,  n  >  R.  One  property  of  (6.2)  can  be  predicted  iitmedlately.  We 
recall  that  (U.23)  represents  simply  a  consolidation  of  (U.lU).  The  present 
sums  J2    ^r>         (reducible  and  irreducible  diagrams  included)  clearly  do  not 
exceed  the  corresponding  sums  in  the  true  problem.  Since  (2.10),  which  consti- 
tutes (U.ll;)  for  that  problem,  is  absolutely  convergent  for  all  t,  it  follows 
that  the  power  series  expansions  of  the  solutions  of  (6.2)  are  absolutely  con- 
vergent for  all  t. 

Nevertheless,  these  solutions  are  not  valid  higher  approximations  to  (2.8). 
The  reason  is  that  the  functions  to  which  their  expansions  converge  become 
infinite  as  t  — >  oo .  None  of  them,  therefore  constitutes  a  uniform  approxi- 
mation, and  for  none  of  them  does  G(co)  exist.  We  shall  illustrate  this  for 
R  ■  2.  If  G(p)  denotes  the  Laplace  transform  of  G(t),  then  (6.2)  for  this  case 
is  equivalent  to 

pG(p)  -  1  -  <b2>  [g(p)]2  >  (b2>  2[a(p)]*^  .  (6.3) 

00 

Let  us  asstune  tentatively  that  G(0)  -  G(t)dt  is  finite.     Then 

^0 


[g(0)]2  -   [l  1  v/=T]/2  (b2>    ,  (6.I4) 

OD 

Which  is  inconsistent  with  the  reality  of  0(t).     Therefore,    I      G(t)dt  cannot  be 

0 
finite.     Further  analysis  readily  shows  that,  for  real  ©,    (6.3)  is  inconsistent 

with  Re  {G(-i<o)}     -  0(co"'^),   a  —t>  0,  if  r  is  any  finite  power.     It  follows  that 

0(t)  grows  faster  than  any  power  of  t  as  t  — >  00. 

The  numerical  solutions  of    (6.2)  for  several  values  of  R  are  compared  with 

(2.8)  in  Fig.  5*     As  R  increases,  it  will  be  noted  that  the  approximations 

increase  in  accuracy  for  small  t  but  diverge  faster  at  large  t.     In  this  respect, 

ou»  present  results  resemble  very  closely  those  of  truncating   (2.10)  after  a 
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finite  number  of  terms.     They  do  not  appesu:  to  represent  a  significant  improve- 
ment over  the  latter. 

The  failure  of  the  present  approximations  has  an  immediate  interpretation. 
Divergence  of  G(t)  as  t  — >  oo   is  inconsistent  with  (2,5).     It  follows  that 
the  values  of  the  irreducible  C's  implied  by  (6.2)  are  not  realizable  by  any 
assigment  of  the  0's  consistent  with   (U.U).     Thus  these  approximations  do  not 
correspond  to  any  dynamic  model  in  our  sense. 

If  we  regard  a  stochastic  assignment  of  the  0    a  ^  a  ^°r  ^  — ^  oo   as  a  dis- 

tiribution  over  the  set  of  index  values  o  and  p,   then  the  C.         are  moments  of 

^n,p 

this  distribution  and  there  are  an  infinite  set  of  realiaability  inequalities 
which  they  must  satisfy.  The  values  Cg.^  -  1,  Cg^j   "  0  (all  higher  irreducible 
diagrams)  for  the  random  coupling  model  correspond  to  complete  statistical  inde- 
pendence of  the  0    „       n  '  V    o     a       in  this  sense.  Nonvanishing  values  for  the 
a,p,a-p    a-p,-p,a 

higher  irreducible  C's  imply  statistical  correlation  among  the  0's. 

The  nature  of  realizability  inequalities  for  simpler  statistical  problems 
suggests  that  when  C  -  -la  wide  choice  of  realizable  nonzero  values  can  be 
given  to  the  higher  irreducible  Cp  .  ,  provided  these  values  are  small  enough. 
Let  us  consider  the  assignment 


C-  ,  -  1,  C,  ,  -  a,  C^    -  0  (all  higher  irreducible  diagrams);      (6.5) 
2;1       hfi  tnjp 


where  a  is  a  real  constant.  For  a  «  1  this  gives  (6.2),  R  -  2.  Instead  of 
(6.3),  we  now  find,  in  the  general  case, 

pG(p)  -  1  -  <b2>  [g(p)]2  *  a  <b2)  ^[kp)]^     ,  (6.6) 
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[g(0)]2  -   [l  -   (1  -  ha)^/2J/2a  (b^)  (6.7) 

Eq,    (6.6)  yields  real,    non-negative  G(a))  for  all  «  if 

-  1/12  <  a  <  1/U  ,  (6.8) 

which  suggests  that  (6.8)  may  represent  the  range  of  realizabilitv  of  (6.5)« 
In  Fig.  6,  the  solution  G(co)  -  n'     Re  { 0(-ico)j  on  the  relevant  branch  of  (6.6) 
is  compared,  for  several  values  of  a,  with  (2.7)  and  (5.6).  The  form  of  G(a)) 
changes  continuously  with  a  up  to  the  limit  a  »  l/h,  where  the  slope  at  co  -  0 
changes  abruptly  from  0  to  oo .  It  is  apparent  that  none  of  the  present  approxi- 
mations represents  a  substantial  improvement  over  the  random  coupling  solution 
(a  =  0).  For  a  <  0,  the  form  of  G(co)  changes  continuously  with  a  down  to  the 
limit  a  -  -  l/l2;  then,  a  singularity  appears  at  the  cut-off  point  of  the  spectnim. 
We  conclude  tentatively,  lacking  contrary  evidence,  that  (6,8)  does  represent  the 
range  of  realizability  of  (6.5). 

The  results  of  the  present  Section  suiggest  that  great  caution  be  exercised 
in  carrying  out  partial  summations  of  diagrams  in  the  power  series  expansion  for 
G(t).  It  is  by  no  means  true  that  the  more  terms  summed,  the  better  the  approx- 
imation. Our  inadmissible  approximation  (6.2)  (R  -  2)  is  equivalent  to  the 
retention,  in  (U.lb)  for  the  true  problem,  of  all  terms  whose  diagrans  can  be 
reduced  to  the  diagram  for  C  ,  by  iterated  application,  on  any  solid  lines,  of 

2J-L 

the  contraction  operations  shown  in  Fig.   7a.     Examples  are  shown  in  Fig.   7b. 
Thus,   the  terms  retained  are  selected  according  to  well-defined  and  plausible 
topological  properties  of  the  diagrams.     Moreover,   as  we  have  noted,   in  the  t 
domain  they  constitute  an  absolutely  convergent  sub-series  of  an  absolutely  con- 
vergent series. 
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It  will  have  been  recognized  by  this  point  that  the  diagram  summations  we 
have  employed  are  intimately  related  to  summations  of  perturbation  series  terms 
in  quantum  field  theory  and  quantum  statistical  mechanics.  Our  present  results 
suggest  that  cairtion  be  exercised  in  these  problems  also.  There  too,  it  is 
possible  that  plausible-appearing  and  summable  classes  of  diagrams  are  better 
omitted  than  included.  We  hope  to  return  to  these  questions  in  a  later  paper. 

7.   SECOND  STOCHASTIC  MODEL 

The  results  of  the  last  Section  ariphasize  the  desirability  of  seeking 
higher  approximations  to  (6.1)  which  correspond  to  realizable  values  of  the  C's. 
We  shall  now  describe  a  second  stochastic  model  for  which  G(co)  satisfies  (2.1i) 
and  is  substantially  closer  to  (2.7)  than  is  the  randcm  coupling  result.  Con- 
sider the  contraction  operation  shown  in  Fig.  8a.  Each  application  to  a  diagram 
reduces  the  nvimber  of  vertices  by  two.  Let  us  take  C-  ^  ■  1  and  assign  the 
value  a"~  to  all  irreducible  C,  ,  whose  diagrams  can  be  transformed  into  that 
for  Cp  ,  by  n-1  applications,  anywhere,  of  this  operation.  These  diagrams 
represent  an  infinite  subset  of  the  terms  in  (I4.23).  (Examples  are  shown  in 

Fig.  8b.)  Let  us  assign  the  value  zero  to  all  other  irreducible  C    .  (Examples 

2n;p 

are  shown  in  Fig.  8c.)  Now  let  us  take  a  =  1.  Clearly  this  implies  the  value  1, 
as  in  the  true  problem,  for  all  reducible  and  irreducible  C„    whose  diagrams 
can  be  transformed  into  that  for  C-  ^  by  repeated  application  of  the  line  opera- 
tion  of  Fig.  Ua  and  the  vertex  operation  of  Fig.  8a.  All  other  C„  ^  have  the 
value  zero. 

We  have  not  found  an  explicit  construction  for  this  model  of  the  type  pro- 
vided by  (5»1),  et  seq.,  for  the  random  coupling  model.  Consequently  we  have 
no  proof  of  realiz ability.  However,  examination  of  the  deperience  of  G(co)  on 
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a,  as  in  Section  6,  suggests  that  the  model  is  realizable,  as  we  shall  see. 
Since  the  present  model  retains  an  infinite  subset  of  terms  in  (U.23), 
it  does  not  directly  yield  a  closed  equation  for  G(t).  However,  we  can  obtain 
a  closed  system  [eqs.  (7.6)  and  (7.1?)]  by  considering  simultaneously  the  first 
two  equations  of  an  hierarchy  analogous  to  (2.1].).  Let 

H(^^)  »Ep<Ha-?,p,a(^)>  '  ^^(tVdt  . 
From  (h.6)  we  find 

dH(t)/dt  -  -  Ep  ^a,p,a-pCp,.p,a  <V-p^a,a(^)>  *  '^^^^        "(°)  "  °  »   ^^.l) 
where 

Since  G^  ^(t)  is  statistically  sharp  (M  — >   00),  it  follows  from  previous  rela- 
tions that  the  first  term  on  the  right  side  of  (7.1)  may  be  rewritten 
-Cg  -j^  (b  )G(t)  .  Hence,  >rtien  Cg.-j^  -  1,  we  have  the  equations 


dG(t)/dt  -  H(t),  0(0)  -  1  , 

dH(t)/dt  =  -(b^)  G(t)  +  J(t),         H(0)  -  0  .  (7,3) 

By  (U.23)  we  have 


00    irr 

n"l   p         '*^ 
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An  analogous  expansion  for  J(t)  may  be  found  by  analysis  very  similar  to  that 
which  leads  to  (U,23).  The  result  is 

n=2   p      "^"'P 

J 
where  ^  is  defined  as  follows:  Construct  a  diagram  part  as  shown  in  Fig.  9. 

Call  it  the  fixed  part.   [The  two  vertices  correspond  to  the  two  0  factors 

J 
which  appear  explicitly  in  (7.2)t]  The  summation  ^  then  is  over  every  closed 

diagram,  constructed  by  combining  two  or  more  vertices  with  the  fixed  part, 

that  contains  no  part,  wholly  external  to  the  fixed  part,  which  is  connected  to 

the  rest  of  the  diagram  bv  only  solid  lines.  The  three  simplest  diagrams  included 

J 
are  shown  in  Fig,  10.  It  will  be  noted  that  ^  contains  both  reducible  and 

irreducible  diagrams  in  the  sense  of  Section  h*     The  reducible  diagrams  all 
resemble  Fig.  10a  in  that  the  associated  C  is  the  product  of  just  two  irreduc- 
ible C's. 

Let  us  denote  by  Fig.  11a  the  totality  of  possible  diagram  parts,  with  the 
two  solid  and  single  dashed  external  lines  shown,  that  can  be  transformed  into 
a  single  vertex  by  repeated  contractions  as  shown  in  Fig.  8a.  Let  us  call  this 
structure  a  consolidated  vertex  part.  Then  we  may  represent  the  entire  class  of 
irreducible  diagrams  which  cwitribute  to  H(t),  in  the  present  model,  by  the 
single  consolidated  diagram  shown  in  Fig.  lib.    In  a  similar  fashion,  we  may 
represent  by  consolidated  diagrams  all  the  contributions  to  J(t)  in  the  present 
model.  Clearly,  two  consolidated  diagrams  included  are  those  shown  in  Figs.  12a 
and  12b.  If  we  independently  replace  the  several  consolidated  vertex  parts  in 
these  diagrams  by  all  possible  diagram  parts  which  they  represent,  then  we  obtain 
all  the  individual  contributing  diagrams  which  are  contract ible  into  Figs.  10a 
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or  10b,  respectively.  The  required  contractions  do  not  involve  the  fixed  part. 
It  can  be  seen  that  all  other  diagrams  contributing  in  the  present  model  are 
represented  by  the  infinite  class  of  consolidated  diagrams  shovm  in  Fig.  12c. 
These  diagrams  all  may  be  transformed  into  that  for  C    by  sequences  of  con- 
tractions  which  now  involve  the  fixed  part. 

As  the  representation  by  consolidated  diagrams  suggests,  the  present  model 
leads  to  a  closed  expression  for  J(t)  in  terms  of  H(t)  and  G(t).  It  is  con- 
venient at  this  point  to  work  with  the  Laplace  transform  representation.  If 

G(p),  H(p),  and  J(p)  denote  the  respective  transforms  of  G(t),  H(t),  and  J(t), 

Ifi 
then  the  transforms  of  (7.3)-(7.5)  are 

pG(p)  »  1  +  H(p)  , 

pH(p)  -  -  (b^)  G(p)  +  J(p)  ,  (7.6) 

H(p)-r''   i:'"'(-l)%„.p<b')"[G(p)]'"  ,  (7.7) 

n«l   p         *^ 

n-2   p        *^ 
Let  U3  write 

J(p)  =  J^^\p)  +  J^2)(p)  ^  3(3)(p)  ^  (^^9) 

where  j'      (p),   J^     (p)»   and  J^^   (p)  are  the  contributions  associated  with 

Figs.  12a,   12b,   and  12c,  respectively.     It  can  be  seen  with  the  aid  of  th«  dia- 

^(l) 
grams  that  the  terms  of   (7.8)  included  in  J^   '(p)  are  in  one-to-one  correspon- 
dence to  the  totality  of  terms  in  the  expansion 

[8(P)]^  -  r        e'"  (-1)"%„,p=2„,;,  <"'>  "*"[8(P)]^"*^    ,  (7.10) 

n,m       p,q 
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which  follows  from  (7.7 )•  Evaluating  the  corresponding  terms  by  means  of  the 
rule  C-    »  a""  (all  nonvanishing  irreducible  C's)  and  then  summing,  we 
readily  find 

J^^^p)  -  [q(p)]-^[h(p)]2  .  (7.11) 

A  similar  correspondence  exists  for  J^   (p),  and  we  thereby  find 

j(^'(p)  .  a[S(p)]-l[H(p)]2  .  (7.12) 

The  evaluation  of  J^  (p)  is  somewhat  more  involved.  Let  us  write 

i-1 

where  J     *      (p)  is  the  contribution  from  the  ith  consolidated  diagram  in  the 
sequence  of  Fig.   12c.     The  series  for  J^-'*    "^(p)  is  in  one-to-one  correspondence 
with  that  for   [h(p)]^,   and  we  find 

j(3.1)(p)  .  ..2  (fc2>  -1[S(p)]-3[h(p)]''  (7.U) 

by  ccnparing  corresponding  tenna.  In  a  similar  fashion,  comparing  the  series 
for  j(3'i*l'(p)  and  [h(p)]2j(3.1)(p),  we  fl«i 

J(3,i*l)(p)  .  .,  (b2>  -1[S(p)]-2[h(p)]2J(3'1'(p)  ,       (7.15) 
Whence, 

J»'(p)  -  j(3.1)(p,  .  ,  (b^)  -l[S(p)]-2[fi(p)]2j(3)(p)  .       (^.is, 


-  37  - 

Combining  our  results  (and  suppressing  the  argument  p)  we  have,  finally, 

Eliminating  H  and  J  from  (7.6)  and  (7.I7),  we  obtain 

(b^>  ^G^  ♦  p(  <b^>  ♦  ap2)G^  -  (  <b^>  ♦  3ap^)G^  ♦  3apO  -  a  -  0  .     (7.18) 

The  relevant  branch  of  (7.18)  is  the  one  for  which  G(co)  •  n"  Re  {G(-iw)}  reduces 
to  (5.6)  when  a  ■  0,  and  it  is  easily  verified  that  such  a  branch  exists.  Let 
us  examine  the  behavior  of  G((a)  on  this  branch  as  a  is  varied.  For  sufficiently 
small  a,  it  is  plausible  to  assume  that  our  assignment  of  values  to  the  C's  is 
realizable.  As  we  increase  a  ,  we  may  plausibly  anticipate  that  G(co)  will 
begin  to  violate  (2.U),  or  at  least  will  exhibit  some  discontinuous  change  in 
its  dependence  on  a,  when  a  critical  value  is  reached  for  which  the  C  »s  are 
unrealizable  by  any  assignment  of  values  to  the  Jjf's  consistwit  with  (U.U).  This 
argument  is  supported  by  the  example  of  Section  6« 

The  behavior  of  (7.18)  is  accessible  by  standard  techniques  for  quartic 
equations.  One  finds  that  the  form  of  G(co)  on  the  branch  of  interest  varies 
continuously  with  a  for  -lA  <  a  <  od  .  Over  this  range, 

G(0)  -  (/?nb^)-^[l  *  (1  *  Ua)l/2jl/2  ^  ^^^^^^ 

and,  for  a  >  0,   G(co)  decreases  monotonically  with  increase  of  00     up  to  a  cut-off 
frequency  given  by 

co/  -     ^a'V^d  *  s)2[l  -  (1  -  Uas)^/2|  ^^2^     ^  ^^^^^ 
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ifhere 

3  -  2  -  (1  .  a-1  .  d1/2)1/3  .  (1  .  a-1  -  0^/2)^/3  , 

D  -  (1  +  a-^)2  -  (1  -  |a-^)3  . 

As  a  increases,  co   increases  monotonically.  For  a  <-l/k,   G(0)  is  complex,  and 
we  conclude  that  this  gives  a  lower  bound  to  the  range  of  realizability.  There 
appears  to  be  no  upper  bound,  at  least  on  the  basis  of  the  present  considera- 
tions. 

The  properties  just  listed  appear  to  justify  the  tentative  conclusion  that 
the  present  model  is  realizable  for  a  -  1.  It  should  be  emphasized,  however, 
that  the  arguments  given  do  not  constitute  a  proof  of  realizability.  The  latter 
would  require  an  explicit  prescription  for  constructing  jf's  >rtiich  yield  the 
model.  We  have  not  found  such  a  prescription,  and  therefore  we  regard  the  inves- 
tigation of  the  present  model  as  incomplete. 

The  function  G(a))  for  a  =  1  is  ccnnpared  in  Fig.  13  with  (5.6)  and  (2.7). 
It  will  be  noted  that  the  present  model  gives  a  close  approximation  to  (2.7) 
and  represents  a  substantial  improvement  over  the  random  coupling  model. 

The  apparent  realizability  of  the  present  model,  and  the  significant  improvo- 
ment  it  gives  over  the  random  coupling  model,  suggest  that  there  may  be  an  infin- 
ite sequence  of  closed,  realizable  stochastic  models  in  lAiich  successively 
broader  classes  of  irreducible  C's  are  given  the  value  1  and  such  that  G(t) 
converges  rapidly  to  its  value  in  the  tnie  problem.  If  so,  the  coupling  coef- 
ficients kr         -1  in  (U»2)  will  exhibit  a  distribution,  as  functions  of  n,  r, 
[n,r,sj 

and  s,  which  clusters  more  and  more  closely  about  the  diagonal  values 

kr-         -,  "  6   5    as  one  ascends  the  sequence.  Thus  we  may  hope  that  any  given 
Ln,r,8j    r,n  s,n 
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dynamical  properties  of  the  models  should  converge  in  a  statistical  sense  to 
those  of  the  original  collection  of  uncoupled  oscillators.     The  analytical  com- 
plexity of  models  higher  than  the  two  already  described  is  formidable. 

8.     DRIVEN  RANDOM  OSCILLATOR 

Let  us  add  to  the  rigltt  side  of  (U.2)  a  forcing  term  f       (t)  \rtiich  ia 
identically  distributed  for  each  n,   statistically  independent  for  different  n, 
and  statistically  independent  of  hr  i  for  all  n  and  r.     Let  us  take  the  initial 
conditions 

^[n](*0^-°>  ^M^^O^-O-  (8-1) 

Writing 

and  introducing  collective  forces  f  (t)  defined  in  correspondence  to  (3.1),  we 
find 

(}r^\it))    -F(t),      (|M-^/2fQ(t)  -f(t)|2>  -0(M-^)  , 

<fJt))-0,       (f^(t)fp*(f)>  -  5^^pF(t,t.)     (a/0)  ,       (8,3) 

where 

F(t,t.).  (f^^j  (t)f(.„j  ••(»■)>  . 

We  see  that  t^it)   plays  a  special  role;  When  M  — >  co ,  the  quantity 
M"^'^fQ(t)  -  M"-"-  YZ    ^f-n]  (*)  ta'<^es  the  sharp  value  f  (t).  It  is  easily  verified 
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that  f   (t)  depends  only  on  the  f .  n  '(t)  if  a  j(  0. 
a  {Jij 

In  place  of   (U.l)  we  now  have 
dq„(t)/dt   .  .iM-1/2  r     9^,,p,^., 


"'/'  ^  "  ^Vp(*)  "  ^a(*>  •  ^^'^^ 


Suppose  that  (8.U)  is  formally  integrated  from  t^  and  solved  by  iteration.  By 
arguments  similar  to  those  which  give  (U.7)  we  find 

<q^(*'>>  -0    (a/0),      (q^(t)q^*(f)>  -0    (a  /  y)  .    (8.^) 
By  (8,5)  and  (3.3)  we  have 

(q^„^(t)>  -q(t),  i(t)  =  (M"^/\(t))  .  (8.6) 

It  may  be  verified  from  a  term-by-term  examination  of  the  iteration  solution 

that  the  variance  of  }f'^^\it)  -  M"^  Y1    ^[n]^*^  ^^  ^^"""""^   *     '^^'^^*  M"^'^^qQ(t) 
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is  statistically  sharp  (M  — *>  oo  )  and  may  be  identified  with  q(t). 

Because  of  the  special  role  played  by  qQ(t),  it  is  convenient  to  impose, 

20 
in  addition  to  (U.U)  and  (U,13),  the  condition" 

^    .       -  1  (jjL  or  cr  -  0)  .  (8.7) 

Then  from  (8.U)  we  obtain  (M  —>  oo  ) 

dq(t)/dt  +  iM"^  r     (b  -q_(t))    -  f(t)  ,  (8.8) 

dq^(t)/dt  -  -ib^q(t)  .  iM-1/2  T    (^,,,.p,pb^.pqp(t)  *  f„(t) 

P 

(a  /  0)  ,  (8.9) 
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where  ^   implies  that  p  ■  0  is  excluded. 
Let  us  write 

q|.^^(t)  -  q(t)  *  qf-nj'(t)  ,  (8.10) 

where,  by  (8.6),  (q,_  »(t)>  -  0.  We  shall  call  q(t)  and  q,  «(t)  the  'coherent' 
and  'incoherent'  amplitudes,  respectively.  An  explicit  solution  for  q(t)  is 
readily  obtained.  From  the  definition  of  the  response  matrix,  and  the  statis- 
tical independence  of  the  f 's  and  the  b's, we  have 

t  t 

a  J  .  *  0  7.* 


^mence,  by  (fi.3),  (8.6),  and  (U.9),  we  obtain 

I 


q(t)  -  I   G(t-s)f(8)ds  .  (8.11) 


Consider  now  the  covariance  of  the  incoherent  amplitude.  Let 

Using  (3.3),  (8.?),  and  (8.6),  we  find 

(Q,-n,„]  (***')>  -m'^E  exp[-i2n(n-m)aA]  (Q^^^(t,f))  ♦  0(M"^)  .     (8.12) 
Therefore,  if  (Q^  ^(t,t»))  is  independent  of  a  (a  /  0),  we  have  (M  — >  oo ) 

Q(t,f)  -  (q^  „(t,f)>     (a  /O)  .  (8.13) 


a.o 
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Let  us  assvune  hereafter  that   (8.13)  holds.     As  we  shall  see  shortly,  this  will 
be  the  case  when  (li.l3)  and   (8.7)  are  satisfied. 
An  important  statistical  property  is 

(bpb_p.   .   .bj,b_^Q<i^^(t,t'))    -  <b^)  .   .   .  <b^>  Q(t,f)  -  0(M';^)  (8.m) 

(a/0,        Ia|    /Ipl   /.   .   ./Inl)   . 

For  the  case  of  uncoupled  oscillators   (all  Jjf's  ■  l)  this  follows  directly  from 
(3.1)  and  the  statistical  independence  of  b       ,  q       (t)  and  b       ,   q  -     (t)  for 

n  /  r.     In  the  general  case,   it  may  be  verified  for  each  term  of  the  iteration 
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expansion  of  the  left  side  of  (B.lU). 

From  (8.9)  we  have 

aQ(t,t«)/dt  -  S(t,t')  ♦  S^(t,t»)  +  Sj,(t,t')  ,  (8.15) 

where 

S(t,t.)  -E  (Sa,a-p,e(*'^')>      («/0)  , 
P 

Sa,a.p,p(*'*'>  -  -^«"'^  <^a,a-p,pVp<lp(*)C(*'>  *  ^8.16) 

S^,(t,t')  -  -i  <b^q^*(t')>  q(t)      (a  /  0)  ,  (8.17) 

Sp(t,f)  -  (q/(t')f^(t)>         (a/0)  .  (8.18) 

Our  notation  anticipates  the  fact  that  the  expressions  given  for  S(t,t»), 
S-(t,t«),  and  S_(t,t')  are  individually  independent  of  a.  It  should  be  noted 

C  F 

that  the  similar  equation  for  dQ(t,t')/9t«  is  redundant  with  (8.15)  because  of 
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the  property 

Q(t,t')  -  Q*(f,t)  ,  (8.39) 

which  follows  from  the  definition  of  Q(t,t  • )  . 

The  quantity  S_,(t,t')  is  readily  evaluated  by  an  argument  similar  to  that 

which  gave  (8.11).  We  have 

t« 
<q>')fjt))  -^J   (Q^^p*(t.-s))(fp*(s)f^(t))d8, 

whence 

t' 
Sp(t,t')  -  f   G*(t«-3)F(t,s)ds  .  (8.20) 

^*0 

The  evaluation  of  (S„  ^  ^  Q(t,t')>  parallels  that  of  (h  „  „  (t))  in 
Section  U.  We  expand  qg(t)  and  q^  (t«)  by  iteration  of  the  integrated  form  of 
(8.9),  leaving  q(t)  (which  is  known)  explicitly  in  the  expansion.  Then  we  aver- 
age and  note  the  regularities  imposed  by  the  sum  rule  for  indices  and  the  sta- 
tistical properties  of  the  b's  and  the  f's.  The  expansion  for  q^  (t » )  involves 
factors  55.  If  we  express  these  as  ^  factors  by  (U.h),  we  are  led,  eventually, 
to  the  irreducible  diagram  e3q)ansion 

,   00    irr 

(Sa,a.p,p(*'^')>-  «   L  .^  E    %p(->-P>P)«2n;p(^'*'>  (^-^l) 

(a,p  /O,   |a|  /  Ipl  /  |a-p|)  , 

where  the  $_       (t,t»)  are  independent  of  o  and  p.     Fran  this  we  obtain 
<in;p 


00         irr 

:  ,E 

n»l      p 


S(t,f)-E    .E        %p«2n,p(*'*')-  (8.22) 
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The  5„   (t,t')  may  be  determined  by  the  variational  method  used  for  the 
2njP 

^   (t)  in  Section  h.  For  variations  which  leave  the  C^    unchanged,  we 
have,  in  correspondence  to  (h.l6), 

^      irr 

The  variation  (U.l?),  with  the  notation  change  p  — >  a-p,  produces  additional 

tenns  on  the  right  sides  of  the  equations  of  motion  for  q  (t ' )  and  qQ('t)«  I" 

-I/2 
correspondence  to  (U,19),  we  find  (to  order  M  '  ) 


t 

Then  (to  order  M"-'")  we  find,  noting  A (2^    o  q  "  '^'S^i 
"  t' 

t 

.j   (Gp^p(t-a)b^.pbp^q/(t.)q,(s)> 


ds  . 


p,p-a,a 


ds 


(8.2li) 


(8.25) 


-  ks 


Using   (8.m)  and  the  sharpness  of  G        and  G        to  reduce  the  averages  in  the 


a,a     ?, 


limit  M  — >  00,  we  have,  finally. 


«2^^(t,t.)-<b2> 


I   G*(t'-3)Q(t,s)d3  -  j   Q(t-s)Q*(t»,8)d3 


(8.26) 


It  is  noteworthy  that  this  expression  depends  on  the  driving  forces  only  impli- 
citly, through  their  effect  on  Q(t,t'). 

The  higher  ^       (t,t")  may  be  found  by  introducing  more  general  variations. 

The  result  is  that  £0    consists  of  a  sum  of  terms  each  of  which  involves  a 
2n;p 

(2n-l)-fold  time  integration  over  a  product  of  2n-l  factors  G,  one  Q  factor, 
and  n  factors  (b  )  • 

We  have  finally  to  evaluate  S-,(t,t')»  It  can  be  shown  from  the  iteration 
solution  of  (8.9)  that  (\<if^it'))   ,  like  (Q^  ^(t,t'))  ,  has  an  irreducible 
diagram  expansion  and  is  independent  of  a  (a  /  0).  From  the  latter  fact  we  hav» 


<b^q/(t.)>  -yr^r^  (Vp*(''^>  *°(«"^)- 

H«nce,  in  the  limit  M  — >  00   we  obtain  from   (8.8)  the  result 
Sc(t,f)    --q(t)[dq*(f)/df   -r*(f)]    . 


(8.27) 


have 


Let  us  now  specialize  to  the  random  coupling  model.  By  (5.3),  we  then 
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S(t,f)  -  42.i(t,f )  . 


(8.28) 


Eqs.  (8.15),  (8.19),  (8,20),  (8.27),  and  (8.28)  now  permit  the  determination  of 
Q(t,t'). 
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Twice  the  real  part  of  (8.15)  for  t  -  t«  is  the  equation  for  the  rate  of 
change  of  the  mean  intensity  Q(t,t)  of  the  incoherent  oscillation.  The  quan- 
tities 2Re|Sp(t,t)]  and  2Re  [Sp(t,t)]  represent  contributions  to  dQ(t,t)/dt 
from  interaction  with  the  coherent  oscillation  and  incoherent  driving  forces 
respectively.  From  (8.26)  and  (8.28)  we  have 

Re  (s(t,t)}  -  0  .  (8.29) 

Thus,  using  (8.27),  we  find 

~  [Q(t,t)  ♦  q(t)q*(t)]  -  0  ,  (8.30) 

when  all  the  driving  forces  vanish.  This  consistency  property,  and  (8.29) 
itself,  are  assured  in  advance  because  our  equations  constitute  the  exact  des- 
cription of  a  model  for  which  (U.b)  holds.  It  is  also  assured  that  the  solution 
Q(t,t»)  of  our  model  equations  obeys  all  the  realizability  conditions  to  which 

22 
covariances  are  subject. 

The  simplest  solution  of  (8.15)  for  the  random  coupling  model  results  from 

taking  forces  which  vanish  except  for  impulses  at  t  -  t^  such  that 

Q(tQ,  tp)  -  1,  q(tQ)  -  0.  Then  it  is  easily  seen  from  (8,26),  and  the  property 

G(t)  -  0*(-t),^^  that  (8.15)  becomes  identical  with  i$,h)  (for  t,t«  >  t^)  under 

the  substitution 

Q(t,t')  -  G(t-f)  .  (8.31) 

Thus  (8.31)  is  the  solution  for  random  shock  excitation,  as  is  required  for 
consistency.  In  general,  Q(t,t»)  and  O(t-t')  do  not  have  the  same  form. 
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9.     NON-QAUSSIAN  FREQUENCY  DISTRIBUTICN 

We  shall  now  descirlbe  briefly  the  generalizations  required  when  the  dis- 
tribution of  br-  -,   is  non-Gaussian.     The  case  (h  r  -y)  /  0  can  be  treated  by 
methods  similar  to  those  used  for  (f  r -.)    /  0  in  Section  8.     But  it   is  simplnr 
to  maintain  the  condition  (b,--|)  -  0  and  instead  replace   (R.I4)  by 


(d/dt  .   ice  .  v)qjt)    *  iM-1/2  Z     ^a,p,a-pNVp(^^    '  ^a^^^    >  (9.1 ) 

P 

where  co  is  the   (real)  mean  frequency  and  we  have  also  included  a  real  damping 
factor  V,      It  is  easy  to  see  that  this  generalization  implies  only  minor  changes 
in  our  treatment  if  b.-  n    remains  Gaussian.     The  factor  t     /2ni  in  the  iteration 
series   (h.m)  is  replaced  by  0^°^  (?K  G^'^b^"  ,   where 

G^^^t)  »  exp[-(iJ^+  v)t]    .  (9.2) 

However,  the  irreducible  diagram  expansions  for  H(t)  and  S(t,t')j  and  the  .expres- 
sions for  the  C,   (t)  and  the  4-  .„(t,t'),  are  unchanged  in  form.  The  effect 
of  w,  V  /  0  is  implicitly  expressed  by  the  changed  values  of  the  functions  G 
and  Q  which  appear  in  these  expressions.  The  only  further  changes  in  Sections 
2-8  are  the  obvious  replacements  d/dt  — >  (d/dt  ♦  i©  +  v)  and  d/dt  — > 
(3/dt  +  ico  +  v)  where  appropriate. 

Now  let  us  consider  the  general  (non-Gaussian)  case  where  the  b  p  -,  are 
identically  distributed  (with  zero  mean)  for  all  n  and  statistically  independent 
for  different  n.  It  is  easy  to  verify  fran  (3.1)  and  (3.2)  that  (3.8)  remains 
valid.  In  place  of  (3.10),  we  find 
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*     ^a^.Y^e«''(<^^)    -3<b'>')       ,  (9.3) 


In  general  we  find   (M  — >  oo  )  that  all  moments  of  the  b     with  indices  equal 

only  °- 

and  opposite  in  pairs  depend  on     (b    )  and  have  the  same  values  as  for  a  Gaussian 

distribution  of  b ^  ,    with  this  variance.     Moments  for  which  the  indices  do  not 

pair    (we  shall  call    them  skew  moments)   have  values  irtiich  depend  on  the  cumulants 

of  the  b  p  -,    distribution;  they  tend  individually  to  zero  as  M  — >  oo. 

The  presence  of  skew  mcments  results  in  new  classes  of  terms  in   (U.ll) 

and  consequently  in   ()i.23).     The  simplest  new  term  in   (U.23)  is 

('±)\.-^   (b^)  G))CG*0,  (9.b) 

where 

%1  '  ""^  IIp^/a,p,a-p^a-0,Y,a.p-/a-p-Y,-p-Y,a  (9.?) 

It  may  be  represented  by  Fig.  lUa.     The  further  terras  represented  by  Figs,   llib 
and  Hic  are  proportional  to     (b^)    -  3  (b  )      and  to     (b  )  (b^)  ,   respectively. 

It  is  possible  to  generalize  otir  sequence  of  models  so  that  closed  equa- 
tions are  produced  vriiich  systematically  include  more  and  more  of  the  information 
expressed  by  the  cumulants  of  the  b  distribution.     We  shall  not  attempt  this 
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here.  However,  it  is  important  to  note  that  the  equations  for  the  random  coup- 
ling model  are  identical  to  those  already  given  no  matter  what  the  (zero-mean) 
b  distribution  may  be.  It  is  clear  that  D^  ^  vanishes  for  this  model  (M  — >  (x>  ), 
and  it  can  be  seen  that  all  the  higher  new  terms  in  (U.23)  vanish  also.  The 
statistical  properties  of  the  random  coupling  model  thus  depend  only  on  the 
variance  (b  )  .  Recalling  (2.3),  which  is  exact  for  the  original  uncoupled 
oscillators,  we  see  that  in  certain  respects  the  random  coupling  model  actually 
will  provide  a  better  approx3.mation  for  distributions  of  b  which  resenble  (5»^) 
in  form  than  for  a  Gaussian  distribution. 

In  the  physical  analogs  to  the  random  oscillator  which  are  our  eventual 
interest,  the  distribution  of  the  stochastic  quantity  corresponding  to  b  may 
itself  be  determined  by  dynamical  processes.   In  this  case  there  may  exist  an 
alternative  to  the  general  treatment  mentioned  above.   It  may  be  physically 
reasonable  to  assume  Gaussian  initial  conditions  for  the  quantities  corresponding 
to  b  and  q.  If  the  dynamical  equations  for  these  quantities  are  then  treated  as 
a  simultaneous  set,  the  non-Gaussian  diagrams  will  not  arise  in  any  of  the  rele- 
vant sequence  of  models.  We  shall  give  an  illustration  at  the  end  of  Section  11, 

10.   PARTICLE  IN  A  RANDOM  POTENTIAL 

Let  the  Schrodinger  equation  for  a  particle  be 

(d/dt  -  iV^)  ^(x^t)  =  -iv(x)  J^(x,t),  (10.1) 

where  v(x)  is  a  real  potential  wMch  is  statistically  distributed  over  an  infi- 


nite ensemble  of  realizations  of  the  system.  This  problem  is  an  exact  homolog 
to  the  randcTO  oscillator  with  respect  to  treatment  by  stochastic  models.  Let 
us  consider  a  collection  of  M  systems  such  that  the  individual  potentials 
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Vr  T  (x)  are  identically  distributed  for  all  n  and  are  statistically  independent 
Ln]  vn 

for  different  n.  Let  ^j.  (x,t)  be  the  Schrbdinger  function  for  the  nth  system* 
Then  we  may  define  the  collective  quantities  (f_(x,t)  and  v  (x)  in  correspondence 
to  (3.1)  and  consider  the  model  equations 


0/6t  -  iV^)  ^  Jx,t)  -  -iM-1/2  ^^  (2(^^p^^_pVp(x)  ^^.p(x,t)  .        (10.2) 

The  ^    a       a   will  be  identical  quantities  for  corresponding  models  in  the  present 
problan  and  the  random  oscillator  problem. 

The  condition  (U.b)  serves  to  maintain  hermitacity  in  the  present  cag-?.  It 
is  easily  verified  from  (U.U)  and  (10.?)  that  the  total  probability 


Ejl^fj.'-'l'd'^'  "  Ejl(»a(l.'-)l'^ 


^'\ 


and  the  total  energy 


/  [■  ^„   hA^^f\-}    *  ^„,.,,  nn]**[n,r.s]    \rf 


['1 


A, 


Where  A-        -i    ^^  defined  by   (U.3),   are  conserved.     Hov;ever,   the  individual 

quantities       [    |  ^p  ,  |   d  x  are  not  constants  of  motion,   in  general.     The  systems 

in  the  collection  exchange  particles  as  well  as  energy. 

Let  us  now  take  the  case  where  v^.  (x)  has  a  multi-variate  Gaussian  distri- 

[nj   «, 

bution.     This  implies  that  all  odd-order  moments  vanish  and  that  all  even-order 

manents  are  expressible  in  terms  of  the  covariance  V(x,x')  •»  (v^  ,  (x)v       (x'))  . 

«v  ^  \    [n]    -     [n]   V.     / 

In  the  collective  representation  we  have 
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(v^(x)vp(x«)v^(x«'  )vjx"'))     -  5^+p\+eV(x,x')V(x'>  ,x"«) 

+  5^^^5p^J(x,x"  )V(x',x'")  +  6^^^6p^V(x,x'")V(x',x"  )  ,  (10.3) 


The  analog  of   (3.fi)  holds,   of  course,   whatever  the  distribution. 

Let  us  define  the  Green's  function  G^      -.  (x,t|x','t')  as  the  solution  (  for 
all  t)  of  the  model  equation  for  if/r  -.  (x,t)  under  the  initial  condition 
^p,(x,t')  «  5       5(x-x'),   and  make  a  corresponding  definition  for  G^     (x,t|x',t'). 
Then,   in  correspondence  to  the  analysis  in  Section  h,   we  find,  when  (U.13)  is 


satisfied, 

00         irr 
"p 


0/3t  -  ±^^)G{x,t\^^,tn  -  n        11        ^-^^\n;p^2n,-p^^**l2^''*'^  ' 


where 


G(x,t'Ix',t')  -  5(x-x')   ,  (lO.li) 

(Q       (x,t|x',t')>     -  6       G(x,t|x',t')   ,  (10.5) 

and  ^^     ^(x,t!x',t')   is  the  hornolog  of  ^„.^(t-t')   , 
cxifP  -^      —  cn,p 

The  functions  J^    .    (x,t|x',t')  may  be  determined  by  the  variational  method 

used  in  Section  )4.     The  variation  ()4.17)  produces,   in  correspondence  to   (U.19), 

the  variation 
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t 
X    [-iM-l/2AC3^^^,v.p(y)0^^^(y,sIx.,t.)J^  (10.6) 


and  we  are  led,  thereby,   to  the  result 
t 


^  ,(x,t|x',f)  «•   f    ds     dV  V(x,y)Q(x,tIy,8)G(y,8!x',f)  , 


(10.7) 


which  corresponds  to  (U.21).  This  result  depends  on  the  fact  that  G   (x,t|x',t') 
13  statistically  sharp  (M  — >  oo  ),  which  ma^  be  demonstrated  in  the  same  way  as 

for  G   (t)  of  Section  U, 

a, a 


The  result  for  ^i^.^(x,t|x',t')  is 


S^^3(x,t|x»,t')  -  j    dsj     dSgj     ds3Jjjd^y3^d3y2d^y37(x,y2)V(y^,y, 


^  °^i^**lll''l^'^^Il*'l'l2'^2^'^fe*^2'j3*°3^^fc'^3'i^'**'^   *  ^^°*^^ 

The  structure  of  expressions    (10.7)  and   (10,8)  may  be  represented,   as  in  Fig.  15, 
by  an  appropriate  labelling  of  the  vertices  in  the  diagrams  for  C„  ,   and  Ci    -c 
The  expiressions  for  all  the  higher  ^   ^    (x,t|x',t')  may  be  written  down  by 
analogy  from  the  diagrams  for  the  corresponding  irreducible  C„^,„. 

CTif  p 

In  the  random  coupling  model,  where  G-  ,  ■  1  and  all  the  higher  irreducible 

C„    vanish,  (10. U)  becomes 
2n;p 


$3 


O/dt  -  iV^2)G(x,t|x',f )  - 


-  J  da  \d\   V(x,y)G(x,t|y,s)G(y,3|x',f)  , 

G(x,t»|x',t')  -  6(x-x')  .  (10.9) 

It  should  be  pointed  out  that  this  result  is  independent  of  the  assumption  that 
the  potential  has  a  Gaussian  distribution,  provided  that  (v  - -,  (x))  •  0  (cf.  Sec- 
tion 9).^"^ 

We  are  assured  that  the  solutions  of  (10,9)  will  exhibit  certain  consistency 
properties  because  this  is  an  exact  equation  for  a  realizable  model.  In  parti- 
cular, if  (10.2)  is  transformed  into  the  raomentuin  representation,  it  follows 
from  a  straightforward  extension  of  the  arguments  given  in  connection  with  (5.7) 
that  G(x,t|x',t')  satisfies  a  basic  spectral  condition.  The  latter  takes  its 
simplest  form  for  the  homogeneous  case  V(x,x')  •  V(x-x'),  in  which  G(x,t|x',t') 
can  depend  only  on  x-x'  and  t-t'.  If  we  write 


Gj^(t-f) 


-  fd^y  G(x,t|x',t')exp(-ik.3r)     (y  -  x-x')  , 

00 

Gj^(co)  -  (2n)"^  I   dsGj^(s)exp(ia)3)  ,  (10.10) 


then  the  spectral  condition  is 

\i<^)   -  |G,,(co)I  .  (10.11) 

jWe  may  note  that  a  (t-t  • )  is  the  diagonal  response  function  for  the  amplitude  in 
the  mode  k.Q  Eq.  (10.11)  implies  the  reciprocity  relation 

G(x,t|x',f)  -  G*(x',f|x,t)  ,  (10.12) 
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also 


which  may  be  deduced  directly  from  (10,2)  and  (U.I4)  in  the  general  inhomogeneous 


When  V(x,x')  =  V(x-x|),  (10.9)  has  the  transform 

t 
0/at  +  ik2)Q^(t)  -  -  f  ds  J  d^k'Vj^_j^,Gj^,(t-s)Gj^(s)  , 


\(0)  -  1  , 


(10.13) 


where 


\  -  (2«)'^  |dV(y)exp(-ik.j;) 


It  is  possible  to  solve  (10.13)  easily  for  very  high  k  (the  WKBJ  limit).  This 
is  of  particular  interest  because  it  is  well  known  that  the  perturbation  approach 
breaks  down  in  this  limit.  Let  us  take  k  sufficiently  high  that  Vj^_j^,  ^  0  unless 
Ik-k'l  «  k.  Then  it  is  plausible  that  we  may  replace  Gj^,(t-3)  by  Gj^(t-3)  in 


(10.13)  and  thereby  obtain 
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where 


t 
(d/dt  +  ik^)Oj^(t)  -  -  (v^)  J  Gj^(t-s)G^(s)ds,     Gj^(O)  -  1  ,     (lO.lU) 
^0  *"     ~ 

(v^)-  V(0)  -  j  Vj^d\. 

Eq.  (lO.lU)  becomes  identical  with  (5»U)  under  the  transformation 


Gj^(t)  -*>  exp(-ik^t)G(t)  , 


<v2>-.><b2>. 


Consequently,   we  have 


\{<^)  -  \ 


(nv^)" 


(ir-'k 


.2\2|l/2 


(|a>-k  I   <  2v^)  (^Q^^^j 

(Ico-k^l   >2v^), 
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where  v^  "  (v  )    .  This  correspondence  between  the  WKBJ  limit  and  the  random 
oscillator  is  not  confined  to  the  random  coupling  model.   If  (10, U)  is  written 
in  the  k   representation  in  this  limit,  and  the  previous  transformation  is  made, 
the  resulting  equation  is  identical  with  (Li,23),  A  particular  consequence  is 
that  the  WKBJ  solution  to  the  true  problem  (all  0's  -  1)  is 


G,.(co) 


(2n(v2>  )-^/2^xp[^l/2(a>-k2)2/  (y^  >J  ,  (10.16) 


in  correspondence  to  (2.7).  This  result  states  that  for  sharp  kinetic  energy 
(sharp  k)  the  total  energy  dlstiribution  follows  the  Gaussian  pxDtential  energy 
distribution.  Considered  as  an  approximation  to  (10.16),  the  random  coupling 
result  (10,15)  exhibits  the  qualitative  physical  fact  that  sharp  momentum  states 
are  not  sharp  energy  states.  The  quantitative  form  of  (10.15)  suggests  that  the 
random  coupling  model  may  represent  a  better  approximation  to  the  original  problem 
if  the  potential  distribution  has  a  clipped  rather  than  a  Gaussian  form  (.cf.  Sec- 
tion 9).  It  should  be  noted  that  the  cumulant-discard  approximation  scheme,  when 
applied  in  the  W?(BJ  limit,  yields  expressions  for  exp(ik  t)G,  (t)  which  are  iden- 
tical in  form  to  (2.13).  This  implies  discrete  spectra  G.  (co),  which  is 
unphysical  compared  to  the  random  coupling  result. 

The  general  correspondence  between  the  WKBJ  limit  and  the  random  oscillator 
includes,  of  course,  the  second  stochastic  model,  discussed  in  Section  7.  The 
WKBJ  results  for  the  two  stochastic  models  and  for  the  true    problem  in  the 
Gaussian  case  are  given  by  Fig.  13,  if  the  horizontal  and  vertical  axes  are 
relabelled  (co-k  )/v^  and  nv^G  (co),  respectively.  Away  from  the  -WKBJ  limit,  the 
analysis  for  the  second  model  is  considerably  more  difficult  than  for  the  random 
oscillator,  although  the  same  in  principle.  The  equations  are  not  reducible  to 
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algebraic  form,  and  the  analogs  to  the  inverses  [B(p)]'  and  (b  )  "  ,  which 
appeared  in  Section  7,  must  be  defined  by  integral  equations. 

Let  us  assume  that  the  Schrodinger  fields  are  switched  on  at  t  «  t-  in  such 
a  way  that  the  W^  |(^j*n^  ^^^   statistically  independent  and  identically  distributed 
for  all  n  and  statistically  independent  of  v   (x)  for  all  r.  Let 

We  shall  caliy/(x,t)  the  coherent  wave  and  [f/  r^^A^j^)   "the  incoherent  wave.  The 
evolution  of  the  coherent  amplitude  and  the  incoherent  covariance  may  be  determined 
by  direct  correspondence  to  the  analysis  in  Section  8. 

Noting  that  our  switch-on  is  equivalent  to  the  action  of  impulsive  sources 
fp  -,(x,t)  "  ^n  -|(x>tQ)5(t-tQ),  we  have,  in  correspondence  to  (8.11), 

i//(x,t)  -  JG(x,t|y,tQ)^(y,tQ)d\  .  (10.18) 

Wh«n  the  potential  is  statistically  homogeneous,  (10.18)  has  the  transform 

V7^(t)-G,(t-tQ)j^/^(tQ),  (10.19) 


(//^(t)  -  (2n)~^  j^(x^t)exp(-ik.x)d^x 


In  this  case  the  various  momentum  modes  of  the  coherent  wave  evolve  independently. 
As  our  WKBJ  limit  results  illustrate,  \(t)   has  a  continuous  spectionn  and  there- 
form  vanishes  as  t  — >  oo  .  Consequently,  the  coherent  wave  eventually  is  ex- 
tinguished by  its  interaction  with  the  random  potential. 
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In  direct  correspondence  to  the  results  obtained  in  Section  8,  we  have 

(^o^2'.**Vy*(2J'»*')>  "  ^a,Y^^^'^*^'**'^     (a/0),   (10.20) 
where  Y(x,t;x*,t')  has  the  symmetry  property 

T(x,tjx',t')  -  Y*(x',t«jx,t)  (10.21) 

and  obeys 

(d/6t  -  iV  ^)T(x,t,-x',t')  -  S(xjt;x',t')  +  S_  (x,t;x',t ' )  ,         (10.22) 

with 

CO    irr 
2^?^**'i^'»*'^  ■  n   H      ^2n;p^2n;p^^*'3S'»"*^'^  (10*23) 

and 

Sp(x,tjx',t')  -  -i^(x,t)0/df  +  iV^^)y7*(x',f)  .  (10.2U) 

There  is  no  terra  corresponding  to  S_(t,t')  because  we  have  not  admitted  sources 


The  functions  $_  (x,t;x',t«),  which  are  homologous  to  the  ^^  ,  (t,t»)  in 
(8.22),  may  be  determined  by  employing  our  variational  procedure  and  noting  the 
statistical  property 

-  V(^,x[)  .  .  .  V(x^,x^)Y(y,t,-y',f)  +  0(M"^) 

(a  /O,   |a|  /  Ipl  /  .  .  .  /  lul)  ,         (10.25) 
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which  corresponds  to   (8.1U).     In  particular,  we  find 

t« 


jd\  V(x,y)    I 
t, 


ds  G  (x',t'|y,s)T(x,tjy,s) 


^2.i(x,t;x',t')  -   Id^y  V(x,y) 

t  ^ 

-    f     ds  G(x,t|y,s)t(x',t'jy,s) 

For  the  random  coupling  model, 

3(x^tjx',t')  -  C2.iU,t;x',f) 


(10.26) 


(10.27) 


and  we  have  a  closed  set  of  equations  which  determine  Y(x^tjx',t')  when  the  ini- 
tial functions  Y(x,tQ^x',tj.)  and  Y  (x,tQ)  are  given.  As  vras  the  case  for 
G(x,t [x',t •),  certain  important  consistency  properties  necessarily  are  exhibited 
by  the  solution  Y(x,tjx^',t ' )  for  any  realizable  model.  In  particular,  we  are 
assured  that  Y(x^t;x',t')  satisfies  all  the  realizability  conditions  to  which  a 
covariance  is  subject.  In  the  homogeneous  case,  where  the  spatial  dependence 
of  Y(x,t;x',t*)  involves  only  x-x',  we  must  have 


Y.  (t,t)  =.  |Y.  (t,t)!  , 


(10.28) 


Yj^(t,t')  -  (2n)"^  j  Y(x^t;x',t')exp(-ik.3r)d^y    (y  -  x-x')  . 

When  t_  — 1>  -co   and  a  statistically  stationary  state  has  been  set  up,   so  that 
Y   (t,t«)  -  Y,  (t-f),   we  have,   further, 


Y^(co)  -   !y^(o))|    , 


(10.29) 
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i^ere 


00 
-00 


7^(co)  -   (2n)"-'   1       T^(t)exp(ia5t)dt. 


Twice  the  real  part  of   (10.22)  for  x»,t'   -  x,t  represents  the  continuity- 
equation  for  the  ensemble  mean  of  the  quantum-mechanical  probability  of  finding 
a  particle.     The  left  side  is  the  quantum-mechanical  equivalent  of  the  substan- 
tial derivative  of  the  mean  probability  density  T(x,t;x,t)  in  the  incoherent 
wave.     The  corresponding  quant: ty  for  the  coherent  wave  is  -2Re  {Sp(x,tjx,t) j 
It  is  clear  from    (10,2?)  and   (10.26)  that  Re  { S(x,t;x,t )]  vanishes.     Consequently, 
the  continuity  equation  simply  states  that  a  particle  enters  the  incoherent  wave 
as  it  leaves  the  coherent  wave. 

The  vanishing  of  Re  (s(x,tjx,t)]  exi^resses  the  fact  that  the  direct  effect 
of  the  potential  on  the  particles  is  to  change  their  momentum  rather  than  their 
position.     To  illustrate  this,    let  us  take  ijl  {xjt)  -  0  and  assume  that  the  poten- 
tial and  the  incoherent  wave  are  statistically  homogeneoixs.     Then  from  the 
Fourier  transforms  of   (10.22),    (10.2?),  and   (10.26)  we  obtain 

t 
df^(t,t)/dt  -  2Re  I     ds  jd\'V^_j^,[Gj^(s-t)Yj^,(t,3) 

*0 

-Oj^,(t-s)Tj^(s,t)]    ,  (10.30) 


after  noting  (10.3  2)  and  (10.21).  The  quantity  T  (t,t)  is  the  mean  probability 
density  for  finding  a  particle  with  momenttun  k.  The  right  side  of  (10.30)  there- 
fore is  the  rate  of  transfer  of  particles  to  this  momentum  for  all  other  monenta 
k^*.  It  is  easily  verified  from  (10.30)  that  |Tj^(t,t )d-^k  is  a  constant  of  motion. 
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Suppose  that  the  fields  have  been  switched  on  at  t^  ■  -oo  in  such  fashion 
that  a  stationary  state  exists  at  time  t.     Using   (10.11)  and  (10.29),   the  right 
side  of   (10.30)  may  be  revfritten  so  that  we  have 

0  =  dTj^(t,t)/dt  -    (    dcod\'Vj^_^,  yco)Y^,(w)  -  Gj^,(co)yco)      .  (10.31) 

-00 


We  note  that  the  right  side  of  (10.31)  is  the  difference  of  two  terms  each 

29 
of  which  is  positive.    The  first  represents  an  input  of  particles  to  mode  k^ 

from  other  modes  k'  and  the  second  represents  an  output  to  these  other  modes. 

If  the  excitation  of  mode  k  only  were  to  be  slowly  increased  by  some  outside 

agency,  it  is  clear  that  the  output  term  would  increase  in  magnitude  while  the 

input  term  would  be  initially  unaffected.  Thus,  the  random  coupling  model  exhibits 

a  plausible  tendency  to  restore  statistical  equilibriuia. 

It  will  be  noted  that  (10,31)  is  satisfied  in  general  if 

T^(co)  -  f(co)Gj^(co)  (10.32) 

where  f (co)  is  a  function  indepereient  of  k.     Now  it  can  be  seen  from  their  defini- 
tions that  G.  (co)  is  proportional  to  the  density  of  eigenstates  of  energy  co  avail- 
able to  a  particle  of  momentum  k,       while  Y  (co)  is  proportional  to  the  occupation 
of  such  states  by  particles  of  this  momentum.     Thus   (10.32)  has  the  usual  form  of 
a  single-particle  equilibrium  distribution  law  if  f  (co)  is  a  fiinction  of  co/e 
(©  ■  tanperature)  appropriate  to  the  statistics  of  the  particle.     In  a  later 
paper,   we  shall  deduce  distribution  laws  of  this  fonn  clirectly  from  a  conciition 
of  statistical  equilibrium  under  small  perturbations  in  the  coupling  among 
systems  in    a     collection,   without  appealing  to  probability  distributions  in  the 
space  of  the  eigenstates  (such  as  the  grand  canonical  distribution). 
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11.     TURBULENCE  DYNAMICS 

The  problem  of  turbulence  dynamics  serves  to  illustrate  the  application  of 
our  methods   bo  equations  of  motion  which  are  nonlinear  in  the  dynamic  variables. 
Tn  order  to  keep  the  formalism  as  simple  as  possible,   we  shall  work  here  with 
the  one-dimensional  scalar  analog  to  the  Naviex-Stokes  equation  proposed  by 
Burgers,         The  treatment  of  the  Navier-Stokes  equation  for  an  inccmipressible 
fluid,   which  we  shall  discuss  briefly,   does  not  differ  in  essentials. 

Burgers'  equation  is 

^  .  V  -4^    u(x,t)  -  .  u(x,t)  ^HCxiil  -  .  i    §>  [u(x,t)]2  .       (11.1) 
at  dx  /  dx  2     dx 

The  function  u(x,t)  may  be  interpreted  as  the  velocity  of  an  infinitely  compres- 
sible fluid,  of  constant  kinematic  viscosity  v,  executing  one-dimensional  motion. 
If  V  ■  0,  the  quantities 


r 


U(x,t)dx,  \  I       [u(x,t)]^dK 


-00  -00 

are  both  constants  of  motion.  We  shall  call  them  'momentum'  and  'energy'  res- 
pectively.  (This  is  not  their  accurate  meaning,  however,  on  the  basis  of  the 
interpretation  just  suggested  for  u(x,t)«^ 

If  an  infinitesimal  forcing  term  5f (x,t),  is  added  to  the  right  side  of  (11.1) 


5u(x,t)  -  J  ds  j  dyG  (x,t|y,s)5f(y,8)  , 
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where  the  infinitesimal  Green's  function  G     (x,t|x',t')  obeys 


^  -  V  -4  )  0^3(x,t|x',f) ^  ru(x,t)Gj.^(x,tIx',f)l 


dt  8x 


G     (x,t'|x',t')  -  6(x-x')    .  (11.2) 

In  correspondence  to  the  procedtire  followed  in  Sections  U,   8,   and  10,   let 
VIS  take  a  collection  of   systems  with  velocity  fields  Ur^-,(x,t),   and  Green's  func- 
tions Gr     -,  (x,t|x',t '),  pass  to  the  collective  representation,   and  consider, 
Ln,mJ 

instead  of   (11.1)  and   (11.2),  model  equations  of  the  form 
dt  3xW      "  "  ~B     -'^»--^'     dx 


a,,/x,t|x.,t.)  -  .M-1/2  rtp.a.^  ^  [«p(->*)^a.p,r(-**l-'>^')]> 


3t  a?y      "*'  ~P  -»K*— K  3x 


0       (x,f|x«,t')  -  6       6(x-x')   .        (ll.U) 

As  before,  the  ^'s  are  independent  of  x  and  t  and  the  same  for  all  ensemble-reali- 
zations of  the  collection. 

We  shall  now  impose  upon  the  5^'s  the  three  conditions 

^a,p,a-p  '  ^a,o-p,p  » 
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The  first  is  a  symmetry  convention.  It  does  not  restrict  the  dynamics.    The 

second  insures  that  (11.3)  preserves  the  property 
u^(x,t)  -  u,*(x,t) 
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and  therefore  the  reality  of  the  Up  -,  (x,t).  The  third  is  identical  vfith  (h.h). 
It  insures  that 

00  CO 


is  a  constant  of  motion  if  v  »  0,     The  property 

,t)dx  -  0 


/OO 

_d 
dt 

-00 


/    "rn](->^ 


follows  from  (11.3)  for  any  values  of  the  (j^'s,    provided  the  u  j. -.  (x,t)  vanish  at 
X  ■  00 .   In  correspondence  to  (8,7)  we  shall  also  require 

iZ(„  .  -  -  1         (n,  X  or  a  -  0)  .  (11.6) 

The  additional  conditions  which  the  ^'s  now  satisfy  imply  only  minor  modi- 
fications in  the  diagraimatic  representation  introduced  in  Section  Ij.  Let  us 
Msociate  with  0     ,   and  !2i   ->  ^  the  vertices  shown  in  Fig.  la  and  lb  Tespec- 
tively.  Then  the  niles  for  associating  diagrams  with  C,  ,  (a,p,a-p)  and  C„ 
are  identical  with  those  given  in  Section  h   if  dashed  lines  are  replaced  by  solid 
lines  everywhere.    We  shall  assume  hereafter  that  (U.13)  is  satisfied. 

Let  us  take 


"  [n]  ^^>*0^  -  "^*'*0^  *  "  [n]  '  ^^'*0^  '  (11-7) 

where  the  initial  values  Uj.-.'  (x^t-)  are  identically  distributed,  with  zero  mean, 
for  each  n  and  statistically  independent  for  different  n.  In  correspondence  to 
(3»8),  it  then  follows  that  the  moments  of  the  u  (x,t-)  vanish  unless  the  sum  of 
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indices  is  zero.  Now  suppose  that  (11.3)  is  solved  by  iteration.  Fi-om  this  pro- 
perty of  the  initial  value  moments,  and  the  combination  rule  for  indices  in 
(11.3),  we  find 

(u^(x,t)up(x',f)  .  .  .)  »  0    (a  +  p  +  .  .  .  /  0).       (11,8) 

Similarly,  the  iteration  solution  of  (ll.h)  yields 

<G^  (x,tlx',t'))  -0      (a/r)  (11.9) 

It  follows  immediately  from  (11.8)  that 

(u^^^(x,t)>  -<M-^\(x,t)> 

for  all  n.  In  correspondence  to  the  similar  result  cited  in  Section  8,  it  can 
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be  shown  that  M  '  u-.(x,t)  is  a  sharp  quantity  (M  — >  00).  Let  us  write 

Uj-  .  (X,t)  -  U(x,t)  +  U|.  ,»  (X,t),     U(x,t)  -  (Uj-^,  (x,t)>  .  (11.10) 

We  shall  call  u(x,t)  and  yp-i'  (x,t)  the  mean  and  fluctuating  fields  respectively. 
Identifying  }f^^^v^{x,t)   with  u(x,t)  In  the  limit  M  — >  co ,  and  noting  (11.6), 
we  may  now  rewrite  (11.3)  and  (ll.h)  in  the  form 


at    dx^ 


]   u(x,t)*u(x,t)  ^^2^  .   1  M-^  ^'  -i  (u  (x,t)u  (x,t)>  ,   (11.11) 
/  Bx  a  3x 

p  .  V  -4  )  "a^^'*^  *  ~  ru(x,t)u  (x,t)l  - 
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G„  Jx,t' |x',t')  -  6^  6(x-x')  .   (11.13) 


where  ^   implies  that  p  ■  0  is  to  be  omitted  and  J^       implies  that  both  p  ■  0 
and  a-p  -  0  are  to  be  omitted.  It  should  be  noted  that  (11.13)  has  the  same 
form  for  a  /  0  and  a  "  0,  Eqs.  (11.11 )  and  (11,12)  are  coupled  equations  which 
determine  the  evolution  of  the  mean  and  fluctuating  fields. 

Now  let  us  assume  that  the  distribution  of  the  initial  values  Up-,'(x,t^) 
is  multivariate  Gaussian.  It  can  then  be  shown  from  the  iteration  solutions  of 
(11.12)  and  (11.13),  using  arguments  similar  to  those  in  Sections  h   and  8,  that 
(^M   (x,t)u_  (x',t'))  is  irdependent  of  a  (a  /  0)  and  that  (  G   (x,t|x',t«))  is 
independent  of  a  (all  a).  Then  it  follows  from  (11.8)  and  (11.9)  that 

(\n]'(^>^)>l'^^''^')>  -^mU^-'^'^'^^'^ 

U(x,t;x',t«)  -  (u^(x,t)u.^(xSf)>     (a/0), 

G(x,t|x»,f)  -  (G^^^(x,tlx«,f)>  .  (II.IU) 

In  correspondence  to  our  previous  results,  G   (x,t|x',t')  is  statistically- 
sharp  (M  — >  00  ),  and  the  covariances  satisfy 
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(u^(x,t)u_^(x',t')ug(y,s)u_p(y',x')   .   .   .>    = 
U(x,t;xSt«  )U(y,s;y»,s')   .   .   .  +  0(M"^) 

(q,P,...  /o,  \a\  4  !pI  /  ...)      (11.15) 

[cf.   (8.lU)Jo     It  follows  from  (11. Ih)  that  U(x,tjx',t ' )  has  the  symmetry  pro- 
perty 

U(x,tjx',t')  »  U(xSt«;x,t)   .  (11.16) 

From   (11.12)  and   (11.13)  we  find  that  U(x,tjx',t«)  and  G(x,t|x',t')  satisfy 
equations  of  the  form 

(_^_v-^  )   U(x,tjx',f)  ■•■  —  ru(x,t)U(x,tjxSt')l    -  S(x,t;x',t')   ,  (11.17) 

Vat       3x^/  3x1-  J 

(±.v  -^1   G(x,tIxSf)   +—  ru(x,t)G(x,t|xSt«)l   ■  H(x,t|xSt'), 

\dt      sx'^y  3x  •-  -^ 

0(x,t'|x',t')  -  6(x-x»),  (11.18) 

where  it  follows  from  the  iteration  solutions,  with  Gaussian  Ur-^-.  (x,tQ)^that 
S(x,tjx',t«)  and  H(x,tlx',t')  have  the  forms 

irr 

S(x,t;x.,t.)-E  E    ^2n:p^2n:p^^'*'^''''^  ^^^-"^^^ 

n   p     »t-   >f 


and 


irr 
H(x,t!xSt.)  -E  Yl        %p^n;p(^'*l^'>*'^  •  ^^^-2°^ 


n   p 
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To  complete  the  set  of  equations,  we  may  rewrite  (11,11)  in  the  form 

^  -  V  -^  )  u(x,t)  +  uU,t)  ^(^Siil  -  -  i  -^  U(x,tix,t)  .  (11.21) 

at     dx"^  /  3x       '^  3x 

Eq.  (11.21)  is  the  balance  equation  for  mean  'momentum'  density  and  (11,1?)  for 
x',t'  "  x,t  is  the  balance  equation  for  the  mean  'energy'  density  in  the  fluc- 
tuating field. 

The  functions  $„   (x,t;x',t«)  and  ^   (x,t|x',t')  may  be  determined  by 
the  variational  procedure  of  Section  U,  using  (11.15)  and  the  statistical  sharp- 
ness of  G^   (x,t!x',t').  The  results  for  $_  ,  (x,t;x',t  •)  and  ^  ,(x,t!x',t') 


f     CO 

;   (x,t;x',t')  -1-^     ds    dy  G(x',t'ly,s)  -^  [U(x,t,-y,3)]^ 

+  -^     ds    dy  G(x,t|y,s)  —  |u(x,t;y,s)U(x',f^y,s)]  ,  (11.22) 

dx    J^         J  dy  ■-  -■  ' 


t      .00 

a 


^.^(x,t|x«,t«)  -  -2    ds  I   dyG(x,t!y,s 


tp    -00 


X  ' 


—  |u(x,t;y,s)G(y,s|x',t')1  ,     (11.23) 
ay  L  J 


In  general  $_  ^  (x,tjx',t')  consists  of  a  sum  of  terms  each  of  which  involves 
a  (2n-l)-fold  space-time  integration  over  a  product  of  2n-l  factors  G  and  n+1 
factors  U.  The  terms  comprising  ^   (x,t|x',t«)  each  involve  a  (2n-l)-fold 
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integration  over  a  product  of  2n  factors  G  and  n  factors  U. 

We  shall  illustrate  the  variational  procedure  in  the  preaenh  case  by  out- 
lining the  analysis  for  ^^  , (x,t;x',t «)•  In  correspondence  to  (8,16),  we  mav 
write,  for  the  present  problem. 

It 

S(x,t,-x',f)  .  J3  (S^,|i,a-p(^**'^'»*'^>  (a/0). 


Then  the  iteration  solution  yields 

00         irr 
"P 


-       00         irr 


(a,.6,a-p  /O)  (|a|   /  |p|  /  |a-p|),  (11.25) 

which  corresponds  to  (8,21),  and  is  the  basis  for  (11.19).  Now  consider  the 
variation  (b.l7).  By  using  (11.5)  several  times,  we  find 

A  of       «  AO^       =■  A(^       -  Ao( 

a-?,-?,a     a-p,a,-p     *"p,-a+p,a     '^p,a,-a,+p 

'"-a,-p,-a+p    '^-a,-a  +  p,-p 
Hence,  recalling  (11.3),  ve   find 

t      00 

Aup(x,t)  .  j  ds  j  dy  Gp,B(x,t|y,3){-M-^/2^(^Q^^^_^^p  -^  K^y'^^^Ki^p^^'^'l], 

tp    -00  ^ 

(11.26) 


to 
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order  M         ,  with  expressions  of  the  same  type  for  Au     „(x,t)  and  Au  _(x',t')» 

These  results  correspond  to   (8.?U).     It  is  important  to  note  that  the  perturba- 

-1/2 
tion  terms  are  0(M     '    )  so  that  the  infinitesimal  response  functiona  may  correctly 

be  used  to  find  the  induced  variations.     Now  we  may  express  A(S     _        -(xjt jx',t ' )) 

to  0(M~   )  in  correspondence  to   (8.?^),   reduce  the  averages  by  using   (11, 15)  and 

the  sharpness  of  the  G       ,   and  appeal  to  the  analog  of   (8.23).     Thereby  we  obtain 

the  result   (11.22). 

The  random  coupling  model  for  the  present  problem  in  nhtainpd  bv  assigning 

the  ^'s  as  in  Section  5,  but  with  the  additional  constraints  (11.5)  and  (11.6). 

It  is  clear  that  these  constraints  do  not  affect  (5.3)  in  the  limit  M  — >  oo  . 

Hence,  we  have 

S(x,t;x',f)  -  42.i(x,tix',t'),   H(x,tlx',t')  -  ^.^(x,t  |x',t  •)         (11.2?) 

for  this  model.  These  relations,  together  with  (Il.l6)-(ll.l8)  and  (11.2l)-(11.23), 
form  a  closed  set  which  determine  u(x,t),  U(x,t;x',t ' ),  and  G(x,t|x',t')  in  terms 
of  the  initial  functions  u(x,tQ)  and  U(x,tQjx«,t-)  . 

The  most  essential  difference  between  the  present  equations  and  the  analogous 
ones  for  the  random  potential  problem  given  in  Section  10  is  that  G(x,t]x',t') 
is  not  independent  of  U(x,t;x',t')  and  u(x,t)  in  the  present  casej  all  three 
quantities  now  must  be  determined  simultaneously,  A  further  consequence  of  the 
non-linearity  is  that  u(x,t)  does  not  have  an  expression  analogous  to  (lO.lB). 
The  Green's  function  G(x,t|x',t')   can  only     describe  the  propagation  of 
infinitesimal  disturbances  6u(x»,t').  In  general,  u(x,t)  /  0  even  if  u(x,tQ)  "  0 
everywhere. 

The  Navier-Stokes  equation  for  the  velocity  u(x,t)  of  an  infinite  incompres- 
sible fliiid  of  kinematic  viscosity  v  may  be  written,  after  elimination  of  the 
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pressure  term,  in  the  tensor  form 


where 


and 


(i  -  "^  ')  "i^2J''^  ■  -  i  'i^nn^^^t'n^^'*^""^!!^^^  *  ^l^-^^) 


i^^v;   ^,,,  V  ^^J^^      i^i,n     a^^axja^ 


V"2f(x)  5  !x-yr-f(y)dV 


=  j!iJ"Z'""^^ 


for  any  f .  We  may  treat  the  incompressible  turbulence  problem  in  direct  analogy 
to  the  foregoing  analysis  by  taking  a  collection  of  flow  systems  with  individual 
velocity  fields  u . '''^-' (x,t )  and  considering  the  model  equation 


(^-0  " 


where  the  u.°(x,t)  are  the  collective  velocity  fields.  Conditions  (11.5)  now 
serve  to  insure  momentum  conservation  and  (when  v  -  0)  energy  conservation. 

The  final  equations  for  the  random  coupling  model  which  result  from  (11.29) 
are  similar  to  those  for  Burgers'  equation,  but  more  conplicated.  In  the  case  of 
homogeneous  turbulence,  they  take  their  simplest  form  when  transformed  to  corres- 
pond to  a  representation  of  the  velocity  field  by  spatial  Fourier  modes.  They 
are  then  identical  with  equations  for  homogeneous  turbulence  derived  previously 
by  a  different  method.    The  earlier  derivation  exploited  the  fact  that  the 
Fourier  amplitudes  of  a  homogeneous  field  have  statistical  properties  which 
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closely  resemble  those  of  the  collective  coordinates  used  in  the  present  paper 
(cf .  Section  3).  Unlike  the  present  approach,  which  involves  no  geometrical 
syirunetiy  restrictions  and  which  may  be  extended  to  bounded  flows,  the  earlier 
treatment  is  valid  only  in  the  homogeneous  case.  An  extensive  discussion  of  the 
energy  dynamics  of  the  random  coupling  model  is  given  in  reference  36. 

We  wish,  finally,  to  give  a  very  brief  discussion  of  turbulent  convection, 
which  will  serve  to  illustrate  a  point  raised  at  the  end  of  Section  9.  Let 
^(x,t)  represent  the  zero-mean  fluctuations  in  the  concentration  of  marked  par- 
ticles carried  along  with  an  incompressible  turbulent  flow  which  obeys  (11.28). 
Then  ^(x,t)  satisfies 


U-^SJ^)    y/(x,t)-.u,(x,t)Ml(i^  , 
Vat      J    r    ^  1  ^     ^^^ 


(11.30) 


where  x  is  the  molecular  diffusivity.  If  y/ j- ,  (x,  t )  and'^^{x,t)   represent,  respec- 
tively, the  individual  and  collective  fields  for  a  collection  of  flows,  the  model 
equation  corresponding  to  (11,30)  is 


—  -  ^v^ ]  ^„(x,t)  -  -M"-'^  r  ^.  o .  ou/(x,t)  •  "-p-^    ,    (11.31) 


dt       /  '  -  — p  -»^»— k^  -  --       8x^ 

where  (11.5)  and  (11.6)  are  satisfied. 

The  random  coupling  model  equations  which  result  from  (11,31),  under  Gaussian 
initial  conditions  of  the  form  we  have  taken  before  are 

—  -  xV,^  +  u.(x,t)  ±-\   G(x,tlx«,t«) 

et    ^   i    ax^y  -  - 

aG(x,tIy,s)  aG(y,s|xSf) 


f  ds  jdVij(x,t;y,! 


ax^       ay. 


G(x,t'Ix',f)  -  5(x-x')  ,  (11.32) 
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dt  ^  ^  '^        dx^    '      ^       -" 


._L[     ds     dV..(x,t;y,s)G(x',f  |y,s) 


aY(x,tjy,s) 


A     f^  o  dT(x',t'; 

+  -£_        ds  dV..(x,tjy,s)G(x,t!y,s) : 

ax,  A  ^J-      -  -      -  3y. 


where 


n,(x,t)  -    (u.M(x,t)>,     U,^(x,tix',f  )  -(u,r"].(x,t)ujM'(x«,f  )> 


T(x,t;x',f )   -     (^^  (x,t)    <fJ^  (x',f)>     , 

and  G(x,t|x',t»)  is  the  mean  diagonal  Green's  function  for  the  concentration 

37 
field«-^'     We  have  assumed  (^p-.(x,t))     =  0,   a  condition  which  is  preserved  by 

the  equations  of  motion,     A  detailed  study  of  the  consequences  of  these  equations 

when  the  velocity  field  is  statistically  homogeneous  has  been  made  by  Roberts, 

who  derives  the  equations  for  this  case  by  methods  related  to  those  of  reference 

36.    Another  special  case  has  been  discussed  by  the  present  author."^^ 

In  accord  with  the  discussion  in  Section  9,  the  random  coupling  equations 

for  turbulent  convection  involve  only  u. (x,t )  and  the  covariance  tensor 

Uj^j(x,t;x«,t»),   regardless  of  the  distribution  of  the  fluctuating  part  of  the 

velocity  field.     Suppose,   now,   we  ask  how  the  higher  statistical  structure  of  the 

velocity  field  can  be  incorporated  in  higher  stochastic  models  for  the  convection 

problem.     If  this  structure  were  known  explicitly,  we  could,   in  principle,   insert 

the  associated  cumulants  in  the  non-Gaussian  terras,   of  the  type   (9.b),  which  con- 
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tribiit*  in  the  higher  models.  An  alternative  procedure  is  to  assume  Gaussian 
initial  conilitions  for  both  the  concentration  field  and  the  fluctuating  velocity- 
field  and  then  treat  (11.29)  and  (11.31)  as  a  simultaneous  set,  making  the  ^^'s 
identical  in  the  two  equations.  The  sequence  of  higher  models  for  this  problem 
would  commence  vdth  that  of  Section  7,  and  the  non-Gaussian  diagrams  would  never 
arise.  The  assumption  of  Gaussian  initial  conditions  often  may  be  physically- 
plausible,  particularly  if  the  flow  has  persisted  long  enough  that  the  higher 
statistical  structure  of  the  velocity  field  is  determined  principally-  by  the 
dynamics  rather  than  by  the  cumulants  of  the  initial  distribution, 
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APPENDIX  A 

The  second  of  relations  (3.12)  may  be  written 

where  we  use  /G   (t))  -  G(t).  Let  the  left  side  of  (A.I)  be  expanded  by 
iteration  of  (U.6).  In  each  term  of  the  expansion  of  G„  „(t)  or  G„  „  (t),  the 
sum  of  indices  of  the  b  factors  is  zero  [cf.  argument  leading  to  (U.?)].  Con- 
sequently, /G   (t))  (^G    (t))  consists  of  terms  of  the  form 

(-i)^-V^^*^^/2  ^  (product  of  (^'s) 


where  there  are  r  factors  b  in  the  first  average  and  s  in  the  second.  For  each 
such  term  there  will  be  a  corresponding  term 

(-i)''"V^^*^^/2  Y2  (product  of  ^-s) 

p,...^t,Y,...<T 

X  (bp...b^b.p.^^^.^b^...b^^.^^^.^)t^Vrisi  (A.3) 

in  the  expansion  of    /g       (t)G         (t)^  ,   where  the  product  of  9^»s  is  identical 
for  given  indices  P,,..,y>..»   •     Let  the  0'3  be  bounded.     Then  the  difference  of 
(A. 2)  and  (A, 3)   is  bounded  in  magnitude  by 

j^-(r*s)/2    |(prod,,ct  of  !2f.s)!         r 

X  t'**Vri3i   .  (A.U) 
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It  now  follows  stranghtforwardly  from  (9.3)  (we  take  the  general  non-Gaussian 

case)  that  (A. It)  is  0(M  )  if  [(product  of  ?f's)  |    is  independent  of  M  (M  ^  oo), 

Similar  analysis  establishes  (8.1h)  and  (11,15)  if  the  iteration  solutions  of 

(8.9)  and  (11.12),  respectively,  are  used  to  express  the  equations  in  terms  of 

the  parameters  and  initial  values,  whose  statistical  properties  are  prescribed. 

In  these  solutions  it  is  convenjent  to  let  the  mean  amplitudes  remain  in  the 

expansions  as  parameters. 

The  significance  of  (3.12)  was  discussed  in  the  text.  Eqs.  (8.lU)  and 
also 
(ll,l5)^may  be  understood  quQlitatlvely  as  consequences  of  the  fact  that  the 

dynamical  behavior  of  a  collective  degree  of  freedom  Is  detemined  (K  — >  co  ) 
by  interaction  with  an  infinite  nunber  of  other  degrees  of  freedom.   The 
dynamical  coupling  with  any  given  few  of  the  other  degrees  of  freedom  is  infini- 
tesimal in  the  limit,  and  this  implies  a  corresponding  weakness  of  statistical 
dependence. 


(a)  (b) 

FIG,    I.  VERTICES    REPRESENTING     <^/i,X,<7    AND  <^*Ji,  X.o- 

FOR     THE    RANDOM     OSCILLATOR 


(a)  (b) 

FIG.  2.  DIAGRAMS    FOR    Cg.,  0  ,  C,,  ,  0,  04^2'' -  ^^D   C4_3() 


FIG.   3.       COMPARISON   OF  G(t)    FOR    THE    TRUE    PROBLEM   AND  THE    RANDOM 
COUPLING     MODEL 


FIG.  4.  (a)  CONTRACTION    OPERATION  FOR   THE  RANDOM   COUPLING 
MODEL 


(b)  TYPICAL    DIAGRAMS    INCLUDED    IN    THE    RANDOM 
COUPLING     MODEL 


FIG.  5.    SOLUTIONS   OF  (62)    COMPARED    WITH    G(t) 
FOR  THE    TRUE    PROBLEM   (R  =  od) 


FIG      6        SOLUTIONS   G(u/)   OF  (6  6)  FOR    SEVERAL   VALUES   OF  a 
COMPARED    WITH    G{w)'  FOR    THE    TRUE    PROBLEM 


ia] 


(b) 

FIG   7.    (a)     CONTRACTION    OPERATIONS    FOR  (6  2)  WITH  R  =  2 
(b)     TYPICAL    DIAGRAMS    INCLUDED    BY  (62)   WITH    R=2 


FIG.  8.  (a)  CONTRACTION    OPERATION    FOR   THE    SECOND  STOCHASTIC    MODEL 
(b)^lC)  TYPICAL   DIAGRAMS   INCLUDED   AND    EXCLUDED  IN  THE  SECOND 
STOCHASTIC    MODEL 


FIXED    DIAGRAM    PART   FOR    J{t) 


)     c 


(a)  (b)  (c) 

FIG      10.        THE    SIMPLEST  DIAGRAMS    CONTRIBUTING    TO    J(t) 


(a) 


FIG    II.     (a)    CONSOLIDATED    VERTEX      PART 

(b)  THE    CONSOLIDATED   DIAGRAM     FOR     H(t) 


TV 


O O 

(b) 


FIG-    12- CONSOLIDATED    DIAGRAMS    CONTRIBUTING    TO    J(t)    IN    THE    SECOND 
STOCHASTIC     MODEL 


"\>N 


13.    COMPARISON   OF    Gtu.)    FOR     THE     RANDOM 
COUPLING    MODEL    AND    THE    SECOND 
STOCHASTIC    MODEL  (DASHED    CURVES) 
AND    FOR    THE    TRUE    PROBLEM 


^O^^' 


(b) 


FIG.    14.      DIAGRAMS    FOR   SIMPLE     NON-GAUSSIAN    CONTRIBUTIONS 
TO     G(t) 


x,t)(  .^ , ^(Vj.s 


(y,s)  (VpS;)^         *         ^(Ys.Ss) 

(a)  (b) 


FIG    15         DIAGRAMS     FOR    C  2;i  (x,  1 1  x',  t' )  AND    U;  3 'i-*  '  i-*' ) 


(q) 


FIG.   16       VERTEX    REPRESENTATIONS    FOR    THE  TURBULENCE    PROBLEM 
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Let  a  device  sample  each  oscillator  in  the  collection  in  turn  at  time  intervals  Tj 
proceeding  in  order  of  increasing  n  and  returning  from  the  Mth  oscillator  to  the 
first  to  repeat  the  cycle  continuously.  At  each  sampling  instant  let  the  device 
produce  a  sharp  pulse,  of  strength  proportional  to  q^  ,  (t).  If  x  is  small  enough 
that  many  cycles  are  completed  before  the  q  r  -,  change  appreciably,  the  spectrum 
of  the  pulse  train  will  approximate  a  line  spectrum  with  frequencies  (i/Vlx    and 
amplitudes  proportional  to  the  q  (t). 
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12»  Reversal  of  the  direction  of  the  arrow  on  a  dashed  line  labelled  \  corres- 
ponds to  the  trivial  notation  change  X  — >  -  X  for  the  svunmed  index  X  in  the 

associated  C_  .„(a,p,a-p)  . 
(in,p 

13.  See  Appendix  A , 

iHe  This  model  is  unrelated  to  the  'random  phase  approximation'  jj).  Pines  and 
D.  Bohm,  Phys.  Rev.  8£,  338  {19^2  f\  ,     We  make  no  assumption  about  the  phases 
of  the  q  themselves. 

1?.  In  general,  S^^  4   (n-l)"""^  . 

16.  The  relations  between  irreducible  and  reducible  C's,  to  which  we  have 
appealed  in  discussing  (6.2),  are  not  affected  by  the  vmrealizability  of  the 
C's.  These  relations  may  be  regarded  here  as  formal  implications  of  the  require* 
ment  that  (U.lli)  agree  with  the  power  series  for  G(t)  obtained  by  the  iteration 
solution  of  (U.23). 

17.  It  is  easy  to  see  that  any  diagram  which  is  transformable  into  the  diagram 
for  C-  -  by  contraction  operations  which  involve  the  fixed  veirtex  may  also  be 
so  transfoimed  by  an  alternative  sequence  of  contractions  which  do  not  involve 
this  vertex.  Thus  the  diagrams  included  in  this  consolidated  diagram  are  exhaus- 
tive. 

18.  The  argument  p  in  Q(p),  etc.,  should  not  be  confused  with  the  index  p  in  . 

^2n;p  • 

19 «     See  Appendix  A. 

20.     Eq.    (8.7)  constitutes  a  very  weak  additional  constraint  on  the  ^'s,  and  we 

anticipate  that,   provided  C^.,   -  1,   it  will  not  affect  the  realizability  of 

any  given  set  of  values  for  the  C^   .     when  M  — ■>  oo  .     This  is  verif iably  so  for 

the  random  coupling  model;     The  sums  ^    C^   ^   (a,p,a-p)  are  xinaffected  (M  — ■>  oo  ) 

for  either  a  «  0  or  a  /  0  . 
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21»  The  svuns  over  (^  products  which  occur  in  the  present  case  involve  ^  rather 

than  the  unrestricted  sunination  by  which  the  C^    are  defined  in  Section  U» 

<in;p 

This  does  not  alter  the  values  of  the  sums  in  the  limit  M  — >  oo ,  however. 
22.  Cvunulant -discard  approximations  of  the  type  discussed  in  Section  2  do  not 
possess  this  consistency  property  and,  in  general,  do  not  yield  (8.29).  The 
cumulant -discard  approach  is  less  satisfactory  for  determining  Q(t,t') 
than  for  determining  G(t).  Cf,  R.H.  Kraichnan,  Phys.  Rev.  107,  IbS?  (1957)  for 
a  discussion  of  an  analogous  situation. 
23*  This  property  of  (5.5)  is  directly  implied  by  (2.)i). 

2U.  However,  their  nujnber  increases  correspondingly  so  that  those  sums  over  indi- 
vidual skew  moments  which  contribute  to  physical  quantities  remain  finite  in  the 
limit. 

25.  We  take  units  such  that  -ft  -  1  and  2m  «  1,  where  m  is  the  particle  mass. 

26.  The  two  problems  may  be  regarded  as  formally  identical  if  b  and  q  are  inter- 
preted as  vectors  in  a  function  space  and  a  correspondence  is  established  between 
d/dt  and  9/^t  -  i V  ^ • 

27.  When  the  latter  condition  is  not  fulfilled,  an  additional  term 

iv(x)G(x,t (x»,t»)  appears  on  the  left  side  of  (lO.U),  where  v(x^  is  the  mean 

'*"  '^  "     then 
potential,  and  V(x,y)  is  defined  in  terms  of  the  zero-mean  part  of  the  potential. 

28.  This  procedure  can  be  justified  a  posteriori. 

29.  V.  necessarily  is  real  and  non-negative. 

30.  We  may  note  that  jd^kGj^(-ie"-^)  -  Q(x^-iQ"-^Ix,0)  is  the  partition  function 
per  unit  volume,  apart  from  normalization. 

31.  J.  M.  Burgers,  Adv.  Appl.  Mech,  1,  171  (19U8)  . 

32.  We  have,  in  fact,  assumed  this  condition  in  writing  (ll.it). 
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33«  The  additional  symmetry  properties  expressed  by  (11.5)  result  in  an  ambi- 
guity in  the  formal  expressions  for  the  C_    given  by  the  rules  in  Section  U, 

^n;p 

There  is,   however,   no  ambiguity  in  value. 

3U.     Only  terms  involving  solely  diagonal  elements  of  the  Green's  function  matrix 
are  of  leading  order  (M  — >  oo  ).     Thus,  for  example,  the  variation  in  u       induced 
by  the  additional  terms  in  the  equation  of  motion  for  U-  does  not  contribute 
in  the  limit, 

3$,    We  use  superscripts  to  label  individual  and  collective  quantities  here  in 
order  to  avoid  possible  confusion  with  tensor  indices. 

36.     R.H.    Kraichnan,    J.  Fluid  Mech.   5,  U97   (1959);  Second  Symposium  on  Naval 
Hydrodynamics,   R.  Cooper,   editor   (National  Academy  of  Sciences-National  Research 
Covmcil,  Washington,   in  press).     The  equations  corresponding  to  the  random  coup- 
ling model  are  called  the   'direct-interaction     approximation' equations  in  these 
papers . 

37«     G(x,t|x',t«)d  X  is  the  probability  that  a  marked  particle  introduced  at 
x',t'  is  in  d-^x  at  x,t. 
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I960.] 
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